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MILITARY  HYDROLOGY 

A  QUASI-CONCEPTUAL  LINEAR  MODEL  FOR  SYNTHESIS  OF  DIRECT  RUNOFF 
WITH  POTENTIAL  APPLICATION  TO  UNGAGED  BASINS 


PART  I :  INTRODUCTION 


Background 

1.  Under  the  Meteorological/Environmental  Plan  for  Action,  Phase  II, 
Headquarters,  US  Army  Corps  of  Engineers  (HQUSACE) ,  has  been  tasked  to  imple¬ 
ment  a  research,  development,  testing,  and  evaluation  program  that  will  pro¬ 
vide  the  US  Army  with  (a)  environmental  effects  information  needed  to  operate 
in  a  realistic  battlefield  environment  and  (b)  the  capability  for  near-real 
time  environmental  effects  assessment  of  military  materiel  and  operations  in 
combat.  In  response,  the  Directorate  for  Research  and  Development,  HQUSACE, 
initiated  the  AirLand  Battlefield  Environment  (ALBE)  Thrust  Program.  Under 
this  new  initiative,  technologies  to  provide  the  field  Army  with  the  opera¬ 
tional  capability  to  perform  and  exploit  battlefield  effects  assessments  for 
tactical  advantage  will  be  developed. 

2.  Military  hydrology,  one  facet  of  the  ALBE  Thrust,  is  a  specialized 
field  of  study  that  deals  with  the  effects  of  surface  and  subsurface  water  on 
the  planning  and  conduct  of  military  operations.  In  1977,  HQUSACE  approved  a 
military  hydrology  research  program;  management  responsibility  was  sub¬ 
sequently  assigned  to  the  Environmental  Laboratory,  US  Army  Engineer  Waterway 
Experiment  Station  (WES),  Vicksburg,  MS. 

3.  The  objective  of  military  hydrology  research  is  to  develop  an 
improved  hydrologic  capability  for  the  Armed  Forces  with  emphasis  on  applica¬ 
tions  in  the  tactical  environment.  To  meet  this  overall  objective,  research 
is  being  conducted  in  four  thrust  areas:  (a)  weather-hydrology  Interactions, 
(b)  state  of  the  ground,  (c)  streamflow,  and  (d)  water  supply. 

4.  Previously  published  military  hydrology  reports  are  listed  on  the 
inside  of  the  back  cover.  This  report  contributes  to  the  ability  to  calculat 
streamflow,  which  is  the  basis  for  developing  improved  flood-forecasting  capa 
bilities  for  use  on  ungaged  watersheds. 


5.  Streamflow  synthesis  from  ungaged  basins  has  long  been  a  subject  of 
scientific  inquiry.  A  survey  of  hydrologic  literature  (Dooge  1976,  Singh 
1978)  suggests  three  fundamental  approaches:  (a)  empirical,  (b)  conceptual, 
and  (c)  physically  based.  The  first  approach  comprises  empirical  relations 
for  determining  some  key  characteristics  of  streamflow  hydrographs,  such  as 
lag  time,  peak  discharge-  time  to  peak,  or  hydrograph  duration.  These  rela¬ 
tions  are  developed  by  standard  curve-fitting  methods  based  on  data  from  gaged 
basins  and  are  then  applied  to  ungaged  basins  with  the  hope  that  they  will 
yield  satisfactory  results.  Although  such  relations  can  be  useful  in  particu¬ 
lar  cases,  this  approach  is,  in  general,  not  scientifically  sound  and  is  often 
discarded  in  favor  of  one  of  the  other  approaches. 

6.  The  second  approach  basically  incorporates  what  are  referred  to  as 
systems  analysis  and  synthesis  techniques  (Dooge  1973,  Nash  and  Foley  1982). 
These  techniques  use  spatially  lumped  parameters,  although  attempts  have  been 
made  to  make  them  quasi-distributed  (Singh  1979).  In  other  words,  they  do  not 
explicitly  take  into  account  spatial  variability  of  rainfall  or  runoff,  even 
though  attempts  have  been  made  to  partly  relax  this  restriction  (Singh  1978) . 
The  major  thrust  has  been  to  develop  the  effective  rainfall-direct  runoff 
relationship.  Effective  rainfall  denotes  that  portion  of  rainfall  which 
becomes  direct  runoff,  whereas  the  remaining  portion  is  denoted  as  abstrac¬ 
tion.  Direct  runoff  is  that  portion  of  streamflow  which  is  composed  of  sur¬ 
face  runoff  and  quick  interflow.  It  is  implicit  here  that  the  volume  of 
direct  runoff  is  equal  to  the  volume  of  effective  rainfall.  A  classic  example 
is  the  unit  hydrograph  approach,  which  is  the  hydrograph  of  direct  runoff  at 
the  outlet  of  a  basin  resulting  from  an  effective  rainfall  of  unit  volume  and 
of  given  duration  occurring  uniformly  in  time  and  space.  It  is  always  associ¬ 
ated  with  the  duration  of  effective  rainfall;  that  is,  as  this  duration 
changes,  so  does  the  unit  hydrograph.  Most  of  these  techniques  therefore 
revolve  around  estimating  the  effective  rainfall,  separating  the  streamflow 
hydrograph,  and  employing  a  spatially  lumped  form  (integrated  over  space)  of 
the  continuity  equation  in  conjunction  with  a  storage-discharge  relation.  The 
effective  rainfall  determination  and  hydrograph  separation  are  somewhat  arbi¬ 
trary,  since  these  are  not  well-defined  concepts  and  are  based  more  on  conve¬ 
nience  than  on  physical  realism. 

7.  On  the  other  hand,  geomorphic  techniques  have  recently  been  advanced 
for  hydrograph  synthesis  (Boyd  1978;  Boyd,  Pilgrim,  and  Cordery  1979; 
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Rodriguez-Iturbe  and  Valdes  1979;  Rodriguez-Iturbe ,  Devoto,  and  Valdes  1979; 
Valdes,  Fiallo,  and  Rodriguez-Iturbe  1979;  Gupta,  Waymire,  and  Wang  1980; 

Wang,  Gupta,  and  Waymire  1981;  Rodriguez-Iturbe  1982).  These  techniques  have 
added  a  new  dimension  to  application  of  geomorphology  to  the  effective 
rainfall-direct  runoff  relationship.  However,  they  remain  to  be  tested  on  a 
wide  variety  of  gaged  basins  and  have  yet  to  be  applied  to  ungaged  basins. 

8.  The  second  approach  is  promising  but  has  shortcomings  that  need  to 
be  properly  addressed.  First,  only  that  portion  of  the  hydrograph  attribut¬ 
able  to  direct  runoff  is  synthesized.  Second,  the  concepts  of  effective 
rainfall  and  direct  runoff  are  not  well  defined.  Third,  the  amount  of  rain¬ 
fall  infiltrating  into  the  ground  is  determined  somewhat  arbitrarily,  although 
there  is  increasing  evidence  to  support  the  view  that  infiltration  is  one  of 

j  the  moot  important  factors  affecting  the  streamflow  hydrograph.  Fourth,  spa¬ 

tial  variability  in  basin  characteristics  affecting  infiltration,  detention 
and  depression  storage,  and  runoff  is  not  accounted  for.  Fifth,  spatial  var¬ 
iability  of  rainfall  cannot  be  handled  analytically  in  a  convenient  manner. 

I 

Sixth,  the  parameters  appearing  in  these  approaches  often  have  little  physical 
significance,  or  they  have  yet  to  be  correlated  dependably  with  physical 
measurements . 

9.  The  third  approach  employs,  in  some  form,  principles  of  mathemati¬ 
cal  physics  which  are  the  laws  of  conservation  of  mass,  momentum,  and  energy 
(Woolhiser  1982) .  The  development  of  techniques  associated  with  this  approach 
has  paralleled,  for  the  most  part,  those  of  the  second  approach;  that  is, 
development  of  the  effective  rainfall-direct  runoff  relationship  has  been  the 
major  thrust.  The  consequence  has  been  twofold:  (a)  the  techniques  have  been 
refined  little  more  than  those  of  the  second  approach  and  (b)  they  have  been 

less  than  practical  working  tools.  Their  extension  to  ungaged  basins  is 

i 

neither  convenient  nor  intuitively  acceptable. 

10.  It  should  be  pointed  out  that  a  few  isolated  attempts  have  been 
made  to  abandon  the  concept  of  effective  rainfall  and  to  consider  infiltration 
and  runoff  simultaneously  during  and  after  the  occurrence  of  a  rainfall 
episode  (Smith  and  Woolhiser  1971;  Rovey  and  Woolhiser  1977;  Singh  1976a, 
1976b;  Singh  and  Agiralioglu  1980,  1981a,  1981b,  1981c;  Sherman  and  Singh 
1976a,  1976b,  1978,  1982).  However,  these  studies  have  been  concerned  princi¬ 
pally  with  overland  flow  and  not  with  other  components  of  streamflow. 
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11.  Although  physically  based  techniques  have  been  successfully  applied 
to  analyses  of  streamflow  hydrographs,  their  application  to  hydrograph  synthe¬ 
sis  for  ungaged  basins  has  yet  to  be  made.  The  reasons  are  manyfold.  First, 
these  approaches  have  not  been  systematically  validated.  Second,  parameters 
such  as  friction  factor  have  been  determined  by  data-fitting  and  not  from 
physical  measurements.  It  is,  therefore  not  clear  if  thesp  parameters  can 
indeed  be  determined  from  commonly  available  measurements  and  have  the  same 
meaning  in  the  context  of  streamflow  synthesis  as  they  are  intended  to  have. 
Third,  a  systematic  data  base  has  not  been  developed  for  obtaining  an  objec¬ 
tive  validation  of  these  techniques  and  their  subsequent  application  to 
ungaged  basins.  Fourth,  dynamic  interactions  with  subsurface  flow  components 
of  streamflow  have  been  evaded.  Fifth,  space-time  quantification  of  friction, 
geometric  complexity,  variability  of  rainfall,  and  variability  in  basin  char¬ 
acteristics  affecting  infiltration  and  runoff  characteristics  have  been  some 
of  the  persistent  problems  yet  to  be  resolved  objectively.  It  is  not  clear 
how  much  detailed  accounting  of  these  factors  is  needed  in  streamflow  syn¬ 
thesis.  No  single  study  can  address  all  of  these  and  related  issues. 

12.  Close  scrutiny  may  suggest  that  a  major  breakthrough  in  streamflow 
synthesis  on  ungaged  basins  is  most  likely  with  a  conceptual  approach.  Recent 
studies  on  application  of  geomorphology  to  basin  hydrology,  cited  previously, 
have  blended  geomorphologic  laws  with  modern  hydrological  systems  analysis  and 
synthesis  techniques.  As  a  result,  they  may  be  on  the  verge  of  providing  a 
unified  framework  for  hydrograph  synthesis  for  ungaged  basins.  This  motivated 
the  use  of  a  ^u^ai-ccncepto^l  apprcr.ch  in  this  studv. 

Objectives 

13.  The  objectives  of  this  study  ar<>  (a)  to  develop  a  quasi-conceptual 
linear  model  for  direct  runoff  hydrograph  synthesis  potentially  applicable  to 
ungaged  basins,  (b)  to  test  this  model  on  gaged  basins,  and  (c)  to  develop  a 
computer  code  for  ready  use  by  field  engineers. 

Scope 

14.  The  theoretical  development  of  the  model  is  presented  in  Part  II  of 
this  report.  Part  III  includes  an  explanation  of  the  model  structure  and  a 
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brief  description  of  each  subroutine.  In  Part  IV  applications  of  the  mode]  on 
five  experimental  agricultural  watersheds  are  discussed.  An  illustrative 
example  showing  the  calculations  required  is  presented  in  Appendix  A.  User 
instructions  are  included  in  Appendix  B  and  a  listing  of  the  program  code  in 
Appendix  C. 


PART  II:  QUASI -CONCEPTUAL  FRAMEWORK 


15.  For  implicit}',  the  assumption  is  made  that  the  transformation  of 
effective  rainfall  to  direct  runoff  is  linear  and  time-invariant.  The  problem 
of  direct  runoff  hydrograph  synthesis  then  reduces  to  determining  the  instan¬ 
taneous  unit  hydrograph  (IUH)  utilizing  basin  morphometry.  The  approach  pre¬ 
sented  here  was  developed  by  Rodriguez-I turbe  and  Valdes  (1979)  and 
generalized  by  Gupta,  Waymire,  and  Wang  (1980). 

16.  Let  it  be  supposed  that  an  instantaneous  burst  of  effective  rain¬ 
fall  having  unit  volume  is  injected  into  the  basin.  This  burst  is  composed  of 
a  large  number  of  particles  n  ,  which  are  noninteracting.  Each  of  these  par¬ 
ticles  will  stay  in  the  basin  for  a  finite  period  of  time  T^,  i  =  1,  2, 

3,...,  n  .  T^  can  be  referred  to  as  holding  time  or  travel  time.  If  it  is 
assumed  that  T^,  1  s  i  <  n,  are  random  variables,  then  these  must  te  inde¬ 
pendently  distributed  by  virtue  of  the  assumption  of  noninteraction  of  parti¬ 
cles.  It  may  be  added  that  the  assumption  of  T  ,  i  =  1 ,  2, . .  . ,  n,  being 
random  is  physically  plausible. 

17.  T^,  1  =  1,  2,...,  n,  depends  on  where  the  ith  particle  lands  on 

the  basin  and  as  a  consequence  the  path  it  takes  to  reach  the  mouth.  The  path 
is  uniquely  determined  by  where  it  lands  in  the  basin.  Obviously,  also 

depends  on  many  other  factors  encountered  along  the  path.  The  paths  available 
for  these  particles  to  follow  are  determined  by  the  basin  geomorphology  in 
general  and  the  channel  network  in  particular. 

18.  Let  the  basin  be  of  W’-order.  Then  the  streams  ,  of  order  i, 

i  =  1,  2,...,  W,  are  available  in  the  basin;  clearly,  S_.  denotes  the  i^ 
order  streams.  In  this  approach,  channel  networks  are  ordered  according  to 
the  Strahler  ordering  scheme  (Smart  1972).  A  particle  goes  through  a  number 
ot  states  determined  by  the  structure  of  the  drainage  network  as  it  travels 
from  its  point  of  landing  to  the  outlet  of  the  basin.  These  states  are  com¬ 
posed  of  overland  regions  and  channels  of  different  orders. 

19.  A  channel  state  of  order  i  is  defined  by  C^,  i  =  1,  2,...,  W, 
as  the  collection  (ensemble)  of  all  the  Strahler  channels  of  that  Strahler 
order.  likewise,  an  overland  region  state  of  order  i  is  defined  by  r^ ,  i 
=  1,  2,...,  W’,  as  the  collection  ot  all  the  regions  draining  directly  into 
the  i*"'1  ordered  channels.  Then  each  particle  will  initially  be  found  in  one 
of  the  overland  states  r^,  i  =  1,  2,  3,...,  W,  and  its  movement  will  be 
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governed  by  the  following  rules  as  a  consequence  of  the  Strahler  ordering 
scheme:  (a)  the  only  possible  transitions  out  of  the  state  r^  are  of  the 

form  r  C^,  1  <  ±  <  W  ;  (b)  the  only  possible  transitions  out  of  the  state 
are  of  the  form  -*■  C  ,  j  >i,  2  ^  j  £  W  +  1,  1  <i<W  ;  and  (c)  there 
is  a  state  ,  defined  as  a  trapping  state.  Transitions  out  of  the  trap¬ 

ping  state  are  impossible. 

20.  These  rules  define  a  collection  S  of  paths,  S  =  { s }  or  seS 

which  a  particle  may  follow  through  to  the  trapping  state,  that  is,  the  outlet 

W-l 

of  the  basin.  For  a  basin  of  order  W  ,  there  are  2  possible  paths.  To 

illustrate,  consider  a  third-order  basin  as  shown  in  Figure  1,  W  =  3  .  The 

path  space  S  =  {s^.s^,  s^,s^}  consists  of  the  following  paths: 

path  s ^ :  r^  ■+  -*•  -*•  -*■ 

path  s_:  r,  -*■  C,  +  C„  -*•  C. 

2  1  i  J  4 

path  s^:  r?  C?  ->  -*■ 

path  s^:  r^  -*•  ■+ 

21.  These  specify  the  spatial  paths  of  a  particle  through  a  geomorphic 
network  of  channels  and  overland  regions.  The  travel  time  of  a  particle  must 
therefore  be  specified  by  the  particular  path  it  takes  to  reach  the  outlet. 

The  travel  time  Tg  is  the  sum  of  the  times  spent  by  the  particle  in  the 
various  states  forming  its  path. 


where  T^  is  the  time  a  particle  spends  in  the  state  x  (x  =  r^  or  for 

some  i)  and  M  is  the  number  of  states.  T  is  assumed  to  be  a  random  vari- 

x 

able.  T^  can  have  an  arbitrary  probability  density  function  (PDF),  and  for 

different  states  x  and  y  ,  T  and  T  can  have  different  PDF's.  How- 

x  y 

ever,  T  and  T  are  assumed  to  be  independent  for  x  *  y  .  The  validity  of 

x  y 

this  assumption  seems  plausible  from  a  physical  standpoint. 

22.  If  T  denotes  the  random  time  that  a  particle  spends  in  the 

B 

basin,  then 


T„  —  I  T 

B  *-  J  s  s 

seS 


(2) 
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where  I  is  the  indicator  function  for  the  path  s  :  that  is,  I  =1  if 
s  s 

the  particle  follows  the  path  s  ,  and  I  =0  otherwise.  The  PDF  of  T  , 

S  B 

denoted  by  f  (t)  ,  is  obtained  as  follows: 
b 

23.  Let  A  .  be  the  ratio  of  the  area  of  r,  to  the  basin  area  A  , 
n  i  w 

and  P^  the  proportion  of  channels  of  oraer  i  merging  into  channels  of 
order  j,j  >i,  2  <  j  <  W  +  1  .  Obviously  P^  =  ^  not 

strictly  true  since  a  basin  of  any  given  order  may  outlet  into  a  stream  sev¬ 
eral  orders  higher.  However,  this  is  convenient  and  does  not  affect  the 
model.  Similarly,  P  =  1  .  Then  for  a  path  seS  of  the  form  s 

=  {Xj,  x2,...,  xk}  where  Xj,  x2>...  xfc  e  {Cj,  C2,...,  Cw;  r^  r2>...,  ry}  . 
The  path  probability  function  is  defined  as 


p(s)  =  Axl 


^xl ,x2’ * ’^xk-1 ,xk 


(3) 


It  should  be  emphasized  that  the  paths  are  all  distinct.  Therefore,  the  prob¬ 
ability  of  T  <  t  is 
B 


P  (T 


<  t)  -  Xj  P(T  <  t)  .  p( 


E  '  7Ts  s 


s) 


-Lfx1  *  Fx2  .....  Fxk(t)  .  p(s). 


(4) 


s  =  {x  ^ ,  x? .  xk  } 


where 

t  =  specific  time 

F  =  cumulative  density  function  of  T 
x  x 

*  =  convolution  operation 

Differentiation  with  respect  to  t  on  both  sides  yields 

£*(t>  'E£k1  *  £x2  fxk 

s  eS 


where  f  u;.  -'as  the  PDF  of  T^  .  Gupta,  Waymire ,  and  Wang  (1980)  have 

established  the  equivalence  of  f  ( t )  and  the  IUH,  h(t)  .  Therefore, 

B 
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•  p(s) 


(6) 


h(t>  =Lf,l  *  f*2  *•••*  £xk 
seS 

where  h(t)  is  the  result  of  an  instantaneous  burst  of  effective  rainfall  of 
unit  volume.  If  the  effective  rainfall  takes  places  continuously  for  some 
time,  then  the  direct  runoff  can  be  determined  by  invoking  the  basin  linear¬ 
ity.  Stated  simply,  the  convolution  integral  can  be  employed  as 

t 

Q(t)  =  f  h  (t  -  t)  I  (t)  dx  (7) 

0 

where 

Q(t)  =  discharge  at  t 
I(t)  =  effective  rainfall 

r  =  variable  of  integration 

Thus,  the  direct  runoff  hydrograph  synthesis  reduces  to  synthesis  of  h(t) 
using  Equation  6. 

24.  In  Equation  6  the  path  probability  function  p(s)  can  be  specified 

completely  from  the  drainage  network  morphometry.  However,  specification  of 

f^  cannot  be  entirely  based  on  physical  considerations.  For  simplicity, 

is  assumed  to  be  exponentially  distributed  with  some  parameter  K  >  0  . 

This  is  consistent  with  the  assumption  of  basin  linearity.  Then  f  ,  *  f  „ 

xl  x2 

*...*  f  k  in  Equation  6  becomes  the  k-fold  convolution  of  independent  but 
nonidentically  distributed  exponential  random  variables.  That  is, 

k 

fxl  *  fx2  *•••*  £xk(t>  -EClk  <«> 

i  =  l 


where  the  coefficients  C  ^  are  given  by  Feller  (1971)  as 

Cik  ■  K*1  Kx2---Kxk-1  [("xl  -  Kxi)  •  (Kxl-1  -  Kxi)"’ 
(Kxl+1  -  Kxl)-"(Kxk  -  Kxi)]'1 


1 L 


in  which  K  ,  *  K.  ,  unless  i  =  k  .  Therefore,  the  IUH  is  given  as 
xi  xk 


k 

h(t)  = 

Cik  6XP  (“Kxit)  '  p(s)’ 

seS  i=l 


(10) 


s  =  ( x j ,  x2,...,  xk) 

25.  To  apply  Equation  10,  the  parameters  K  must  be  determined. 
Following  Gupta,  Waymire,  and  Wang  (1980),  the  mean  holding  time  of  an  i 
order  Strahler  channel  (state)  is  given  by 

J-.  a(z.y/\  i  <  i  on 

where  a  is  an  empirical  constant  and  is  the  average  channel  length  of 

order  i  ,  which  can  also  be  computed  as 


L.  = 


z 

j=l 


31 


i  =  1,  2 . W 


(12) 


where  N.  is  the  number  of  streams  of  order  i  and  L . .  is 

th  1  Jl 

the  j  stream  of  order  i  .  Likewise,  the  mean  holding  time 


,  th 


order  overland  region  can  be  given  by 


the  length  of 

1/K  of  an 
ri 


1 


/ArlA.\I/3 


,  1  <  i  <  W 


(13) 


From  a  physical  point  of  view.  Equations  11  and  13  state  that  the  mean  holding 
time  of  a  given  state  is  proportional  to  some  "characteristic  length"  of  the 
state.  The  constant  a  is  determined  empirically  and  plausibly  may  remain 
more  or  less  constant  from  one  state  to  another  within  a  given  basin.  Addi¬ 
tional  work  will  provide  its  range  of  variability  on  basins  of  diverse  geomor- 
phologic  characteristics. 
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26.  To  use  Equations  11  and  13,  the  constant  a  must  be  specified. 
The  first  moment  of  the  IUH,  h(t)  ,  being  equal  to  the  mean  holding  time  of 
the  basin,  K  ,  can  be  written  as 

D 


t  h(t)  dt 


00 

/ 

00 

r 

£  t  Q(t)  dt 

J  t  I(t)  dt 

oo 

r 

00 

r 

/  Q(t)  dt 

/  Kt)  dt 

J 

_0 

J 

_0 

(14) 


(15) 


From  Equations  10  and  14  it  can  be  shown  that 


h  =X  p(s) 

seS 


(16) 


s 


<x, 


x2>  x3,.. 


xk} 


27.  If  Equations  11  and  13  are  substituted  into  Equation  16,  the  only 
unknown  is  a  .  However,  K  is  estimated  following  Boyd  (1978)  as 

D 


*8 


=  b 


.0.38 


(17) 


where  is  in  hours  and  is  in  square  kilometres.  The  parameter  b 

must  be  determined  empirically.  Thus,  for  a  specified  value  of  K  ,  a  can 
be  determined.  Methods  for  obtaining  b  are  discussed  later. 
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PART  III:  A  QUASI- CONCEPTUAL  LINEAR  MODEL 


28.  The  quasi-conceptual  model  based  on  drainage  basin  morphometry  for 
direct  runoff  hydrograph  synthesis  (GMHS)  consists  of  a  number  of  subroutines, 
each  describing  a  unique  component.  The  arrangement  of  components,  as  shown 
in  Figure  2,  depends  upon  the  need  for  optimization  of  model  parameters.  If 
optimization  is  not  required,  the  components  are:  (1)  MAIN,  (2)  BASIN, 

(3)  LAG,  (4)  HOLD,  (5)  IUH,  (6)  PRECIP,  (7)  NEWTON,  (8)  INFIL,  (9)  XDATA,  and 
(10)  CONVOL.  On  the  other  hand,  when  optimization  of  parameters  is  required, 
components  are  (1)  MAIN,  (2)  BASIN,  (3)  EXOP,  (4)  PRECIP,  (5)  NEWTON, 

(6)  INFIL,  (7)  XDATA,  (8)  BROSEN,  (9)  OBJECT,  (10)  LAG,  (11)  HOLD,  (12)  IUH, 
and  (13)  CONVOL.  A  flowchart  of  the  model  is  given  in  Figure  3,  and  its  com¬ 
puter  listing  is  provided  in  Appendix  C.  A  brief  discussion  of  the  subrou¬ 
tines  follows. 

29.  The  component  MAIN  outputs  general  information  on  the  GMHS  model, 
initializes  parameters,  reads  in  and  outputs  the  model  objective,  and  speci¬ 
fies  some  inputs  required  by  subroutines  later.  It  also  monitors  whether 
optimization  of  model  parameters  is  required.  Put  succinctly,  MAIN  sets  the 
stage  for  the  model  and  the  tasks  to  be  performed  by  it. 

30.  The  rainfall-runoff  data  are  processed  by  the  subroutine  PFF.CIP. 
These  data  are  properly  arranged,  and  their  units  specified.  First,  the  rain¬ 
fall  data,  which  include  values  of  rainfall  intensity  versus  time,  are  read. 
Since  time  is  read  in  clock-hours,  it  is  reduced  to  a  time  series.  Runoff 
data,  which  include  values  of  discharge  versus  time,  are  then  read.  Here 
also,  the  time  values  are  reduced  to  a  time  series.  The  runoff  data  represent 
direct  runoff.  If  hydrograph  separation  needs  to  be  performed  for  computation 
of  the  direct  runoff,  a  separate  subroutine  must  be  provided  for  this  purpose. 

31.  Effective  rainfall  and  the  portion  of  rainfall  not  contributing  to 
direct  runoff  are  computed  by  using  subroutines  INFIL  and  XDATA.  The  effec¬ 
tive  rainfall  data  are  properly  arranged.  The  time  difference  between  the 
start  of  the  effective  rainfall  and  that  of  the  direct  runoff  is  noted. 
Infiltration  capacity  is  computed  as  a  function  of  time  using  the  Philip  two- 
term  infiltration  model  (Philip  1969)  .  If  the  infiltration  capacity  is  to  be 
computed  by  another  method,  INFIL  must  be  modified  accordingly.  The  infiltra¬ 
tion  model  has  two  parameters:  (a)  sorptivity  accounting  for  capillary 
effects  and  (b)  saturated  hydraulic  conductivity  accounting  for  gravity 
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effects.  These  parameters  are  computed  in  an  iterative  manner  based  on  New¬ 
ton's  method  with  the  subroutine  NEWTON.  It  is  assumed  that  sorptivity  is 
subject  to  change  from  one  rainfall-runoff  event  to  another  on  the  same  basin; 
on  the  other  hand,  saturated  hydraulic  conductivity  remains  fixed  for  a  basin 
but  may  differ  from  one  basin  to  another. 

32.  The  basin  characteristics  are  analyzed  by  the  subroutine  BASIN. 

The  principal  geomorphologic  characteristics  are  (a)  basin  area,  (b)  areas  of 
overland  regions,  (c)  channel  lengths,  and  (d)  number  of  channels  of  each 
order.  This  subroutine  is  used  to  calculate  mean  channel  lengths  and 

areas  of  overland  regions  for  each  order.  Basin  lag  is  computed  using  basin 
area  in  association  with  Equation  17  by  the  subroutine  LAG.  If  a  different 
method  is  to  be  used  for  computing  basin  lag,  this  subroutine  must  be  modified 
accordingly. 

33.  The  mean  holding  times  of  overland  flow  and  channel  flow  are  com¬ 
puted  by  the  subroutine  HOLD,  using  Equations  11  and  13  and  the  basin  char¬ 
acteristics  given  by  the  subroutine  BASIN.  The  instantaneous  unit  hydrograph 
is  computed  by  the  subroutine  IUH  using  Equations  9  and  10.  To  obtain  the 
direct  runoff  hydrograph,  the  IUH  is  then  convoluted  with  the  effective  rain¬ 
fall  obtained  from  the  subroutine  PRECIP  by  the  subroutine  CONVOL,  which  also 
compares  computed  direct  runoff  hydrographs  with  the  corresponding  observed 
direct  runoff  hydrographs. 

34.  When  optimization  of  parameters  is  needed,  then  some  additional 
components  are  used  as  shown  in  Figure  2.  The  subroutine  EXOP  provides 
pertinent  information  required  by  the  optimization  algorithm,  including 
specification  of  initial  guesses,  upper  and  lower  bounds  on  parameter  values, 
number  of  stage  searches,  and  convergence  limit. 

35.  The  subroutine  OBJECT  specifies  the  objective  function  to  be  used 
in  optimization  of  model  parameters.  The  objective  function  was  defined  as 
the  sum  of  squares  of  deviations  between  observed  and  computed  discharge  peaks 
and  their  times  of  occurrences.  A  weighting  factor  was  used  to  assign  rela¬ 
tive  weights  to  the  two  components  of  the  objective  function.  Not  more  than 
20-percent  weight  was  allocated  to  the  component  based  on  the  sum  of  squares 
of  deviations  between  observed  and  computed  peak  times. 

36.  Optimization  of  parameters  is  performed  by  the  subroutine  BROSEN, 
which  combines  the  original  Rosenbrock  method  (Rosenbrock  1960),  the  Palmer 
version  (Palmer  1969),  and  the  penalty  function  method.  The  problem  of 
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optimization  is  formulated  as  a  constrained  minimization  problem  requiring  the 
vector  always  to  be  an  interior  point  of  the  feasible  set.  The  subroutines 
EXOP  and  OBJECT  provide  pertinent  information  to  initiate  optimization. 
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PART  IV:  APPLICATION  TO  NATURAL  WATERSHEDS 


37.  The  quasi-conceptual  linear  model  presented  previously  was  verified 

on  five  small  experimental  agricultural  watersheds  designated  as  C,  D,  G,  Y, 

2 

and  2-H.  These  watersheds  range  in  area  from  0.0137  to  17.72  km  .  The  avail¬ 
ability  of  rainfall,  runoff,  and  geomorphic  data  was  the  primary  consideration 
for  their  selection.  The  geographic  locations  of  these  watersheds  are  shown 
in  Figure  4.  Watersheds  C,  D,  and  G  are  shown  in  Figure  5,  while  watersheds  Y 
and  2-H  are  shown  in  Figures  6  and  7,  respectively. 

Description  of  Watersheds 


Watershed  C 

38.  Watershed  C  is  located  near  Riesel,  TX.  As  shown  in  Figure  8,  it 

2 

is  a  second-order  watershed  having  an  area  of  2.343  km  .  Its  tree-structure 
is  shown  in  Figure  9.  Its  drainage  network  properties  are  abstracted  from  the 
topographic  map.  The  order  of  channel  network,  number  of  channel  elements  of 
each  order,  and  length  and  area  of  each  channel  element  are  given  in  Table  1. 
Watershed  D 

39.  Watershed  D,  shown  in  Figure  10,  includes  watershed  C.  Located 

2 

near  Riesel,  TX,  it  has  an  area  of  4.492  km  .  It  is  a  second-order  watershed 
having  a  tree-structure  as  shown  in  Figure  11.  Its  drainage  network  proper¬ 
ties  are  shown  in  Table  2. 

Watershed  G 

40.  Watershed  G,  located  near  Riesel,  TX,  includes  watersheds  C  and  D. 

2 

It  has  a  total  area  of  17.72  km  ,  as  shown  in  Figure  12.  This  is  a 
fourth-order  watershed,  as  shown  in  Figure  13.  Its  drainage  network  proper¬ 
ties  are  given  in  Table  3. 

Watershed  Y 

41.  Watershed  Y,  shown  in  Figure  14,  is  located  near  Riesel,  TX,  and 

2 

has  an  area  of  1.251  km  .  This  is  a  third-order  watershed,  as  shown  in  Fig¬ 
ure  15.  Its  drainage  network  properties  are  presented  in  Table  4. 

Watershed  2-H 

42.  Located  near  Hastings,  NE,  Watershed  2-H  is  the  smallest  of  all 

watersheds  considered  in  this  study.  As  shown  in  Figure  16,  it  has  an  area  of 
2 

0.0137  km  .  It  consists  of  three  channel  elements  as  shown  in  Figure  17. 
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Table  1 

Drainage  Network  Properties  of  Watershed  C,  Riesel,  TX 
(Watershed  Area  =  2.343  km2) 


Serial 

Channel 

Length 

Contributing  Area 

9 

Number 

km 

ft 

km 

acres 

Order  1 

1 

1.295 

4,250 

0.833 

205.91 

2 

0.647 

2,125 

0.232 

57.40 

3 

0.610 

2,000 

0.272 

67.19 

4 

0.687 

2,255 

0.257 

63.46 

5 

0.555 

1,820 

0.230 

56.82 

Order  2 

1 

0.88? 

2,895 

0.519 

128.22 
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Table  2 

Drainage  Network  Properties  of  Watershed  D,  Riesel,  TX 
(Watershed  Area  =  4.492  km^) 


Order  1 

1 

1.295 

4,250 

0.833 

205.91 

2 

0.647 

2,125 

0.232 

57.40 

3 

0.609 

2,000 

0.272 

67.19 

4 

0.687 

2,255 

0.257 

63.46 

5 

0.554 

1,820 

0.230 

56.82 

6 

1.143 

3,750 

0.426 

105.40 

7 

1.256 

4,120 

0.450 

111.25 

8 

0.838 

2,750 

0.624 

154.32 

9 

0.480 

1,575 

0.132 

32.51 

Order  2 

1 

2.108 

4,940 

1.006 

168.66 

1,975 

80.11 
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Table  3 


Drainage  Network 

Properties  of  Watershed 

G,  Riesel ,  TX 

(Watershed  Area  =  17.72  km2) 

Serial 

Channel  Length 

Contributing 

Area 

Number 

km 

ft 

km 

acres 

Order  1 

1 

0.765 

2,510 

0.646 

159.72 

2 

1.753 

5,750 

0.989 

244.40 

3 

1.4478 

4,750 

0.620 

153.30 

4 

2.118 

6,950 

1.083 

267.64 

5 

0.363 

1,190 

0.189 

46.92 

6 

0.399 

1,310 

0.101 

24.96 

7 

1.0866 

3,565 

0.757 

187.05 

8 

1.256 

4,120 

0.450 

111.25 

9 

1.143 

3,750 

0.426 

105.40 

10 

0.555 

1,820 

0.230 

56.82 

11 

0.687 

2,255 

0.257 

63.46 

12 

1.295 

4,250 

0.833 

205.91 

13 

0.648 

2,125 

0.232 

57.40 

14 

0.610 

2,000 

0.272 

67.19 

15 

0.838 

2,750 

0.624 

154.32 

16 

0.480 

1,575 

0.132 

32.51 

17 

0.777 

2,550 

0.474 

117.21 

18 

0.686 

2,250 

0.779 

192.53 

19 

0.533 

1,750 

0.233 

57.51 

20 

0.3429 

1,125 

0.097 

23.95 

21 

1.067 

3,500 

0.625 

154.42 

22 

0.839 

2,755 

39.15 

96.73 

23 

0.570 

1,870 

0.310 

76.52 

24 

0.419 

1,375 

0.116 

28.72 

25 

0.968 

3,175 

0.528 

130.51 

26 

1.343 

4,405 

0.680 

168.00 

Order  2 

1 

0.917 

3,010 

0.334 

82.53 

2 

3.216 

10,550 

1.396 

345.06 

3 

0.326 

1,070 

0.077 

19.12 

4 

0.954 

3,130 

0.512 

126.60 

5 

1.646 

5,400 

0.698 

172.53 

Order  3 

1 

3.394 

11,136 

2.521 

623.02 

2 

0.155 

510 

0.077 

19.14 

Order  4 

1 

0.107 

350 

0.029 

7.31 
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Table  4 


Drainage  Network  Properties  of  Watershed  Y 

,  Riesel ,  TX 

(Watershed  Area  =  1.251  km2) 

Serial 

Number 

1 

km 

Channel  Length 

ft 

Contributing  Area 
km^  acres 

1 

0.395 

Order  1 

1,300 

0.282 

69.655 

2 

0.097 

1,450 

0.097 

23.938 

3 

0.332 

1,090 

0.152 

37.673 

4 

0.094 

310 

0.094 

23.349 

5 

0.137 

450 

0.122 

30.020 

1 

0.296 

Order  2 

970 

0.112 

27.606 

2 

0.543 

1,780 

0.216 

50.819 

1 

0.259 

Order  3 

850 

0.2 

41.211 

This  is  a 

second  order  watershed.  Its  drainage  network  properties  are 

given 

in  Table  5. 


Table  5 

Drainage  Network  Properties  of  Watershed  2-H,  Hastings, 
NE  (Watershed  Area  =  0.0137  km2) 


Serial 

Channel 

Length 

Contributing 

Area 

Number 

km 

ft 

km 

acres 

Order  1 

1 

0.0219 

72 

5.79  x  i0'3 

1 .4298 

2 

0.015 

4 

0.001 

0.2468 

Order  2 

1 

0.062 

204 

0.007 

1.7217 
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Rainfall-Runoff  Data 


43.  Rainfall-runoff  data  for  each  watershed  were  obtained  from  the 
US  Department  of  Agriculture  publications  entitled,  "Hydrologic  Data  for 
Experimental  Agricultural  Watersheds  in  the  United  States."  These  publica¬ 
tions  contain  the  largest  yearly  flood  events,  between  8  and  10  for  each 
watershed.  These  events  were  divided  into  two  mutually  exclusive  groups,  one 
for  optimization  of  model  parameters  and  the  other  for  model  verif ication. 
Numbers  of  events  available  for  each  basin  and  used  for  model  calibration  and 
verification  are  as  follows: 


Number 

of  Rainfall-Runoff 

Events 

Available 

Used 

Used 

Watershed 

for  Analysis 

for  Calibration 

for 

Verification 

C 

9 

5 

4 

D 

8 

4 

4 

G 

8 

4 

4 

Y 

8 

4 

4 

2-H 

10 

5 

5 

For  each  rainfall-runoff  event,  direct  runoff  was  obtained  by  hydrograph 
separation . 


Determination  of  Infiltration 


44.  Infiltration  for  each  rainfall-runoff  event  was  determined  on  each 
watershed  by  using  the  Philip  two-term  infiltration  model  (Philip  1 969) , 

f  =  A  +  0.5  St'0,5  (18) 


where 

f  =  rate  of  infiltration  (cm/hr)  at  time  t 

A  =  parameter  approximately  equal  to  saturated  hydraulic  conductivity 
(cm/hr) 

S  =  parameter  called  sorptivity  (cm/hr^’^) 

The  parameter  A  depends  mainly  on  the  soil  type  and  was  therefore  fixed  for 
a  given  basin.  Values  used  were  as  follows- 
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Watershed  Value  of  A,  cm/hr 


C  0.254 

D  0.254 

G  0.254 

Y  0.254 

2-H  0.508 

The  parameter  S  depends  on  antecedent  soil  moisture  and  other  physical 
characteristics.  It  was  determined  for  each  rainfall-runoff  event  on  each 
basin  by  a  volume  balance  analysis.  Its  determination  on  an  ungaged  basin 
remains  an  unsolved  problem. 


Parameter  Estimation 


45.  The  GMHS  has  only  one  unknown  parameter  in  Equation  17.  This 
parameter  b  was  determined  for  each  basin  by  using  a  modified  Rosenbrock- 
Palmer  optimization  algorithm  (Rosenbrock  1960,  Palmer  1969).  The  values  for 
the  various  basins  were  as  follows: 


Watershed 

C 

D 

G 

Y 

2-H 


Value  of  b,  cm/hr 
0.875 
0.875 
1.2734 
0.875 
0.875 


Instantaneous  Unit  Hydrograph 


46.  Using  these  parameter  values,  the  IUH  was  determined  for  each 
watershed;  the  IUHs  are  shown  in  Figures  18-22.  It  is  apparent  that  the  IUHs 
possess  appropriate  shape  characteristics.  For  very  small  watersheds,  C  for 
example  (Figure  18),  the  IUH  experiences  a  quick  rise  and  a  quick  recession. 

As  the  area  increases,  the  rates  of  rise  and  recession  become  more  moderate  as 
can  be  observed  for  watershed  G  (Figure  20) . 
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Runoff  Prediction 


47.  The  runoff  hydrograph  was  predicted  for  each  event  in  the  predic¬ 
tion  set  using  the  parameter  b  ,  estimated  in  the  manner  set  forth  in  para¬ 
graph  44.  Comparisons  of  observed  and  predicted  runoff  hydrographs  for  sample 
events  on  each  watershed  are  shown  in  Figures  23-27.  The  predicted  hydro¬ 
graphs  compare  reasonably  well  with  observed  hydrographs  with  regard  to  shape, 
time  of  rise,  time  of  recession,  and  peak  characteristics.  The  prediction 
error  in  peak  discharge  and  time  to  peak  is  as  high  as  50  percent;  in  most 
cases,  though,  it  is  considerably  less.  Two  factors  are  worthy  of  note  here. 
First,  antecedent  moisture  conditions  are  extremely  important.  The  infiltra¬ 
tion  parameter  S  and  the  effective  rainfall  pattern  are  very  sensitive  to 
the  antecedent  moisture  condition  and,  as  a  consequence,  so  is  the  runoff 
hydrograph.  A  small  change  in  the  effective  rainfall  pattern  results  in  a 
marked  difference  in  runoff  hydrograph  characteristics.  Second,  the  parameter 
b  ,  although  determined  optimally,  may  not  have  represented  the  range  of  con¬ 
ditions  persisting  on  a  given  watershed  over  a  long  period  of  time.  This  is 
due  to  a  relatively  small  number  of  ev>  its  being  available  for  its  estimation. 
The  runoff  hydrograph  is  quite  sensitive  to  b  since  this  is  the  only  param¬ 
eter  in  the  TUH.  Nevertheless,  given  model  simplicity  and  its  basis  in 
drainage  network  morphometry,  the  prediction  results  are  encouraging.  Addi¬ 
tional  model  testing  needs  to  be  done  for  more  definitive  conclusions. 

Considerations  of  Basin  Size 


48.  Although  the  GMHS  has  been  applied  to  five  small  gaged  basins,  its 
application  is  by  no  means  confined  to  small  basins.  Large  basins  have 
pronounced  variability  in  rainfall  distribution,  infiltration  rate,  and  surfi- 
cial  characteristics,  all  of  which  need  to  be  accounted  for  in  the  model. 

There  are  two  ways  to  handle  this  problem.  First,  the  entire  basin  may  be 
considered  one  unit,  regardless  of  how  heterogeneous  it  is.  The  basin  is 
represented  by  a  number  of  paths,  each  having  an  associated  area  corresponding 
to  an  ensemble  of  the  portions  of  basin  area  draining  into  this  path.  Because 
these  portions  are  of  a  heterogeneous  nature,  hydrologic  variables  can  be 
averaged.  For  example,  a  certain  path  drains  some  of  the  overland  regions  of 
first  order.  Rain  falling  on  these  regions  can  be  proportioned  by  their 


respective  areas,  and  the  same  can  be  done  for  infiltration  and  other 
variables . 

49.  Second,  a  large  basin  can  be  divided  such  that  each  subbasin  can  be 
considered  a  homogeneous  unit.  The  model  can  then  be  applied  to  each  sub¬ 
basin,  and  outputs  of  the  subbasins  properly  routed  to  produce  the  direct  run¬ 
off  hydrograph  of  the  entire  basin.  Therefore,  the  size  of  the  basin  does  not 
appear  to  be  a  limitation  on  model  applicability. 
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PART  V:  CONCLUSIONS 


50.  The  following  conclusions  are  drawn  from  this  study: 

a.  The  IUH  determined  by  the  model  appears  to  possess  appropriate 
hydrologic  properties.  From  those  generated  and  examined  in 
this  study,  it  is  apparent  that  they  possess  appropriate  shape 
characteristics.  For  very  small  watersheds,  the  IUH 
experiences  a  quick  rise  and  a  quick  recession.  As  the  area 
increases,  however,  the  rates  of  rise  and  recession  become  more 
moderate . 

b.  The  runoff  hydrographs  predicted  by  the  model  compare  reason¬ 
ably  well  with  observed  hydrographs  with  reference  to  shape, 
time  of  rise,  time  of  recession,  and  peak  characteristics.  The 
prediction  error  in  peak  discharge  and  time  to  peak  was  as  high 
as  50  percent;  in  most  cases,  though,  it  was  considerably  less 
than  30  percent. 

c:.  Antecedent  soil  moisture  and  infiltration  are  extremely  impor¬ 
tant  for  accurate  model  predictions.  The  infiltration  param¬ 
eter  S  and  the  effective  rainfall  pattern  are  very  sensitive 
to  antecedent  moisture  conditions  and,  as  a  consequence,  so  is 
the  runoff  hydrograph.  A  small  change  in  the  effective  rain¬ 
fall  pattern  makes  a  material  difference  in  the  characteristics 
of  a  predicted  runoff  hydrograph.  The  runoff  hydrograph  is 
quite  sensitive  to  b  since  this  is  the  only  parameter  in  the 
IUH.  The  parameter  b  ,  although  determined  optimally  in  this 
study,  was  probably  not  representative  of  the  range  of  condi¬ 
tions  that  persisted  on  a  given  watershed  over  a  long  period  of 
time  because  of  a  small  number  of  events  available  for  its 
estimation. 

d.  The  b  parameter  appearing  in  the  lag-area  relation.  Equa¬ 
tion  17,  needs  further  scrutiny.  This  parameter  should  be 
related  to  some  physical  basin  characteristic. 

e.  The  GMHS  model  is  partially  based  on  drainage  network  proper¬ 
ties.  This  feature  suggests  that  the  model  should  be  appli¬ 
cable  to  ungaged  basins.  However,  additional  model  testing 
will  be  needed  to  make  more  definitive  inferences. 
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PART  VI :  RECOMMENDATIONS 


51.  This  report  represents  a  portion  of  a  larger  effort,  i.e.  the 
simulation  of  streamflow  for  ungaged  basins.  Much  additional  work  is  needed. 
Some  fruitful  areas  of  research  are  as  follows: 

a.  Determination  of  the  volume  of  direct  runoff  resulting  from  a 
specified  rainfall  event  is  essential  for  subsequent  synthesis 
of  its  associated  direct  runoff  hydrograph.  Current  procedures 
for  computing  this  volume  are  inadequate  and  usually  are  not 
applicable  to  ungaged  basins.  Despite  its  crucial  importance 
in  streamflow  simulation,  this  aspect  has  not  been  addressed 
adequately  in  hydrologic  literature. 

b.  A  study  to  relate  the  b  parameter  in  Equation  17  to 
measurable  basin  characteristics  is  required.  This  is  essen¬ 
tial  if  geomorphologic  approaches  are  to  be  used  to  synthesize 
the  IUH  for  ungaged  basins. 

£.  Evaluating  the  effect  of  basin  size  and  its  ordering  on  the  IUH 
is  important  from  a  practical  standpoint.  The  detail  required 
for  describing  a  drainage  network  should  be  determined  for  the 
model  reported  here.  For  example,  is  it  necessary  to  represent 
a  sixth-order  basin  as  it  is,  or  will  scaling  down  to  fourth- 
order  representation  suffice? 

d.  The  effects  of  spatial  distribution  of  rainfall  on  generation 
of  direct  runoff  are  not  completely  known.  This  is  an 
important  aspect  of  streamflow  forecasting  and  deserves  con¬ 
siderable  attention. 

e.  The  sensitivity  of  the  GHMS  model  to  various  kinds  of  errors  in 
its  parameters  and  inputs  needs  to  be  determined.  This  is 
necessary  to  decide  whether  the  model  is  adequate,  requires 
improvement,  or  can  be  further  simplified  without  significant 
loss  of  accuracy. 

f_.  For  the  model  to  be  applicable  to  ungaged  basins,  each  of  its 
components  needs  to  be  related  to  measurable  basin  characteris¬ 
tics.  Parameters  of  the  infiltration  model  might  be  estimated 
in  this  manner. 

g.  A  better  assessment  of  the  accuracy  and  reliability  of  this 
model  is  needed.  The  level  of  confidence  that  can  be  placed  on 
model  results  is  not  clear. 

h.  The  GHMS  model  should  be  compared  with  others  on  the  basis  of 
drainage  network  characteristics.  Results  of  such  an  effort 
will  allow  for  placing  the  model  in  its  proper  perspective, 
especially  in  relation  to  others. 

i.  Based  on  applications  made  to  date,  i.e.  to  small  basins,  the 
GHMS  model  is  best  interpreted  mathematically  in  terms  of  the 
standard  hydrologic  concept  of  storage  elements.  In  the 
future,  though,  when  applied  at  the  subbasin  level  where 
routing  becomes  an  integral  part  of  the  overall  procedure,  the 
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GHMS  should  be  interpreted  mathematically  as  representing  a 
network  of  storage  elements  and  channels. 
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INTERIOR  DRAINAGE  BOUNDARY 


Figure  L.  A  hypothetical  third-order  watershed  with  Strahler  ordering  system 


Figure  2.  Components  of  the  hydrograph  synthesis  model  GMHS 


Figure  3.  Computer  flowchart  of  the  GMHS  model 


Figure  4.  Location  of  watersheds 
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Figure  9.  Tree  structure  of  watershed  C 
Riesel,  TX 
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Figure  11.  Tree  structure  of  watershed  D 
Riesel ,  TX 


Figure  13.  Tree  structure  of  watershed  G,  Riesel ,  TX 
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Figure  15.  Tree  structure  of  watershed  Y,  Riesel,  TX 
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Figure  20.  The  instantaneous  unit  hydrograph  of  watershed  G, 

Riesel,  TX 


Figure  21.  The  instantaneous  unit  hydrograph  of  watershed  Y, 

Riesel,  TX 
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Figure  24.  Comparison  of  observed  and  predicted  runoff  hydrographs 
for  the  rainfall-runoff  event  of  16  July  1961  on  watershed  D, 

Riesel,  TX 
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Figure  25.  Comparison  of  observed  and  predicted  runoff  hydrographs 
for  the  rainfall-runoff  event  of  29  March  1965  on  watershed  G, 

Riesel,  TX 


RUNOFF,  CM/HR 


TIME,  MIN 

Figure  26.  Comparison  of  observed  and  predicted  runoff  hydrographs 
for  the  rainfall-runoff  event  of  24  April  1957  on  watershed  Y, 

Riesel ,  TX 
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Figure  27.  Comparison  of  observed  and  predicted  runoff  hydrographs  for 
the  rainfall-runoff  event  of  3  July  1959  on  watershed  2-K,  Hastings,  NE 


APPENDIX  A:  AN  ILLUSTRATIVE  EXAMPLE 


Watershed  2-H,  located  near  Hastings,  NE,  illustrates  the  highlights  of 
the  quasi-conceptual  model.  This  watershed  has  been  discussed  previously  in 
the  main  text.  The  steps  involved  in  using  this  model  are  given  below. 

1.  Compute  the  watershed  area.  Watershed  2-H  has  an  area  of 
Aw  =  0.0137  km2  . 

2.  Order  the  channel  links  according  to  the  Strahler  ordering  scheme  as 
shown  in  Figure  16.*  Draw  the  subwatershed  boundaries  for  each  channel  link. 
Watershed  2-H  is  a  second-order  basin.  Its  drainage  basin  properties  are 
given  in  Table  5. 

3.  Measure  the  length  and  area  of  each  link  and  overland  region.  For 
watershed  2-H,  these  are  shown  in  Table  5. 

4.  Compute  the  average  values  of  length  and  area  for  channels  and  over¬ 
land  regions  respectively  for  each  order.  For  watershed  2-H, 

Lj  =  0.0184  km 
L  =  0.062  km 
A  =  0.00340  km2 
A  *  0.0070  km2 

5.  Determine  the  path  space  and  the  paths  for  the  watershed.  In  the 

(2-1) 

present  case,  the  number  of  paths  is  2  =2 


by 


and 


i  2 

are  defined  as 


The  path  space  then  is  S  =  {s^ . 


Let  the  paths  be  denoted 
s^).  The  individual  paths 


’  1  ‘ 


T*  C  **  C 

1  1  2 


^  2  *  r2  c2 

6.  Compute  for  each  path 
the  watershed  2-H, 

0.0068 


and 


the  ratio 


Ari.  i  =  1,  2 


For 


rl 


0.0137 

0.007 


=  0.496 


=  0.511 


r2  0.0137 
7.  Compute  the  quantity 
P 


ci,cj 


In  the  present  case, 


c  1  ,c2 


f  =  1 


*  See  figures  and  tables  in  the  main  text. 
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p 


c2,c3 


1 

1 


c 


3 


represents  the  trapping  state. 

8.  Compute  the  path  probabilities  p(S) 

p(s.)  =  A  P  .  P  .  .  =  0.496  x  1.0  x 

r  1  rl  rl,cl  cl,c2 


p(s.)  =  A  P  =  0.511  x  i.o  =  0.511 

2  r2  r2’C2 


These  will  be: 
1.0  =  0.496 


9.  Compute  the  basin  lag.  Tf  b  in  Equation  17*  is  assumed  to  be 
0.875  and  the  exponent  is  0.38,  then 

=  0. 875(0. 0137)0'38  =  0.171  hr 

10.  Compute  the  mean  holding  time  of  each  overland  flow  region 
and  each  channel  order  by  using  Equations  11-13  in  conjunction  with  Equa¬ 
tion  16.  For  the  watershed  2-H, 


=  a  x  0.4520 


1/3 

1  /0.511  *  0.0137  \ 

Kr2  a\2  x  i  x  0.062  ) 

=  a  x  0.3836 


~  =  a(0.0184)1/3  =  a  x  0.2640 
cl 


=  a(0. 062) 1/3  =  a  x  0.3958 
c2 


*  See  equations  in  the  main  text. 
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Therefore 


0.171  =  a[ (0.4520  +  0.2640  +  0 . 3958) (0. 496) 
+  (0.3836  +  0.3958) (0.511) ] 

=  a(0. 5514  +  0.3983) 

=  a(0.950) 


This  yields 


0.171 

0.950 


0. 180 


Using  this  value  of  a  ,  then 


~  =  0.08136  ; 

K.  , 
rl 

~  =  0.06900 
r2 

—  =  0.04752  ; 

cl 

~  =  0.07124  ; 


K  ,  =  12.291 
rl 

K  .  =  14.483 

r2 

K  .  =  21.044 
cl 

K  „  =  14.036 
c/. 


11.  For  each  path,  arrange  values  of  the  inverse  of  the  mean  holding 
time  in  a  vector  according  to  the  elements  involved  in  the  path.  For  the 
watershed  2-H, 


path  s^: 


path  s^: 


<  r 


r  r 


<  K 


rl* 


K  „>  -► 
cz 


-  12.291,  21.044,  14.036  > 

<  r„,  c„  >  ♦  <  K  _,  K  _  >  -►  <14.483,  14.306  > 
L  L  XL  CL 


The  path  probability  vector  is 


p(s)  : 


<  Sj,  S£  >  <  0.496,  0.511  > 
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12.  Compute  the  values  of  C^  for  each  path  and  s2  using  Equa¬ 
tion  9.  In  the  present  case,  the  following  is  obtained  for  the  path  : 


K  ,  K  . 
rl  cl 


12.291  x  21.044 


'13  (Kcl  -  Krl)(Kc2  -  Krl)  (21.044  -  12.291) (14.036  -  12.291) 


=  16.934 


K  ,  K  . 
rl  cl 


12.291  x  21.044 


'23  (K  -  K  ) (K  _  -  K  )  (12.291  -  21.044) (14.036  -  21.044) 

rl  ci  cz  cl 


K  .  K  . 
rl  cl 


=  4.2166 


12.291  x  21.044 


'33  (K  ,  -  K  _)(K  ,  -  K  )  (12.291  -  14 . 036) (21 . 044  -  14.036) 

r  i  cz  cl  cz 


=  -21.151 


and  for  path  s^  , 


r2 


14.483 


'12  (K  „  -  K  _)  (14.036  -  14.483) 

cz  r2 


K 


r2 


14.483 


'22  (Kr2  -  Kc2)  (14.483  -  14.036) 


=  -32.400 


=  32.400 


13.  Compute  the  IUH  using  Equation  10.  For  the  watershed  2-H,  the 
following  is  obtained, 

h(t)  =  [C13  exp  (-K^t)  +  C23  exp  (-K^t)  +  c33  exP  (-K^t)  ] 

p(sL)  +  [C12  exp  (-K^t)  +  C22  exp  (-K^t)]  p(s2) 

=  [16.934  exp  (-12.291t)  +  4.2166  exp  (-21.0440 

-  21.151  exp  (-14.0360]  0.496  +  [-32.400  exp  (-14.4830 
+  32.400  exp  (-14.0360]  0.511 


For  different  values  of  time,  the  instantaneous  unit  hydrograph  (IUH)  can  be 
computed  as  shown  in  Table  Al. 
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Table  A1 

IUH  for  Watershed  2-H  Located  near  Hastings,  NE 


Time 

Time 

h(t) 

min 

hr 

1/hr 

0 

0.000 

0.0000 

2 

0.033 

0.1954 

4 

0.067 

0.2909 

6 

0.100 

0.3124 

8 

0.133 

0.2906 

10 

0.167 

0.2489 

12 

0.200 

0.2021 

14 

0.233 

0.1579 

16 

0.267 

0.1200 

18 

0.300 

0.0893 

20 

0.333 

0.0653 

22 

0.367 

0.0471 

24 

0.400 

0.0336 

26 

0.433 

0.0238 

28 

0.467 

0.0167 

30 

0.500 

0.0116 

32 

0.533 

0.0081 

34 

0.567 

0.0056 

36 

0.600 

0.0038 

38 

0.633 

0.0026 

40 

0.667 

0.0018 

42 

0.700 

0.0012 

44 

0.733 

0.0008 

46 

0.767 

0.0006 

48 

0.800 

0.0004 

50 

0.833 

0.0003 

APPENDIX  B:  USER  INSTRUCTIONS 


1.  The  quasi-conceptual  model  based  on  drainage  basin  morphometry  for 
direct  runoff  hydrograph  synthesis  (GMHS)  requires  data  only  on  storm  rain¬ 
fall,  soil  infiltration  characteristics,  and  the  drainage  network  characteris¬ 
tics  of  a  basin.  Thus,  the  model  can  potentially  be  applied  to  synthesize 
direct  runoff  hydrographs  on  ungaged  basins.  To  obtain  data  on  drainage  net¬ 
work  properties,  it  is  sufficient  to  have  a  topographic  map,  preferably  with  a 
scale  of  1:24,000.  Topographic  maps  for  most  of  the  basins  in  the  United 
States  are  available  from  the  US  Geological  Survey.  Data  on  rainfall  and  soil 
infiltration  characteristics  used  in  this  study  were  obtained  from  the 

US  Department  of  Agriculture  publication  entitled  "Hydrologic  Data  on  Experi¬ 
mental  Agricultural  Watersheds  in  the  United  States." 

2.  The  GMHS  contains  a  number  of  subroutines,  the  use  of  which  depends 
upon  whether  parameter  optimization  is  or  is  not  required.  The  arrangement  or 
sequencing  of  the  subroutines  is  shown  in  Figure  2.*  A  computer  program  was 
developed  and  is  available  in  the  form  of  a  Fortran  IV  deck.  The  major  func¬ 
tions  of  the  program  are  shown  in  Figure  3. 

3.  As  for  all  programs,  the  preparation  of  input  data  is  critical. 

Some  common  requirements  are  as  follows.  All  integer  numbers  must  be  right 
justified,  that  is,  placed  as  far  to  the  right  in  the  available  field  as  pos¬ 
sible.  Decimal  points  are  necessary  unless  integer  numbers  are  used.  When  a 
decimal  point  is  used,  it  must  occupy  a  location  in  the  field  just  as  an 
integer  would.  For  example,  the  number  19.8934  would  require  at  least  seven 
spaces  in  the  field.  If  more  than  one  card  of  the  same  format  is  included  in 
the  deck,  the  location  of  the  decimal  points  should  be  kept  the  same  from  one 
card  to  another  to  facilitate  key-punching  of  the  cards.  The  following  dis¬ 
cussion  provides  information  on  input  variables,  data,  and  formats  for 
specific  subroutines  in  the  program. 


*  See  main  text  for  figures  and  tables. 
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GMRS :  MAIN 


4.  This  constitutes  the  main  program.  It  provides  general  information 
about  the  model,  for  example  its  purpose,  and  calls  for  execution  of  the 
model.  Its  input  is  given  as  follows: 

a.  Specify  the  purpose  of  the  computer  program  such  as,  "The  pur¬ 
pose  of  this  program  is  to  synthesize  a  runoff  hydrograph  using 
drainage  network  properties."  This  statement  is  denoted  by  PURP 
and  appears  on  cards  1  and  2  at  the  start  of  the  program.  The 
user  may  enter  any  alphanumeric  information  on  columns  1-80  of 
two  consecutive  cards.  This  is  specified  as  (PURP(I) ,1=1 ,40) 
using  an  A-format  as  FORMAT (20A4) .  This  information  will  be 
printed  at  the  beginning  of  the  computer  output  to  indicate  the 
purpose  of  the  program. 

b.  Specify  the  time  interval  of  computation  and  the  number  of 
basins  under  study.  These  are  denoted  respectively  by  DT  and 
NW  and  are  given  on  card  3.  The  format  for  reading  them  is 
FORMAT (FI 0.4, 15) . 

£.  Specify  the  number  of  rainfall-runoff  events  for  which  the 

program  is  to  be  used.  This  is  denoted  by  NOBS  and  specified  on 
card  4.  The  format  for  reading  it  is  FORMAT (15). 

d.  Specify  the  parameter  A  of  the  Philip  two-term  model.  This  is 
denoted  by  AA  and  specified  on  columns  1-10  of  card  5.  This 
parameter  is  assumed  constant  for  a  given  basin  but  may  vary 
from  one  basin  to  another.  On  the  same  card  are  specified  EX 
and  NXM,  which  denote  the  exponent  of  the  lag-area  relationship 
and  the  number  of  time  intervals  of  computation.  These  are 
entered  into  columns  11-20  and  21-25,  respectively.  The  format 
for  reading  all  three  of  them  is  F0RMAT(2F10.4,I5) . 

e.  A  control  designated  as  KOPT  is  given  on  card  6.  An  integer 
number,  either  0  or  1 ,  is  specified  and  determines  whether 
optimization  of  model  parameters  is  or  is  not  required.  When 
KOPT  is  0,  optimization  is  not  needed.  When  it  is  1,  optimiza¬ 
tion  is  needed.  KOPT  is  entered  into  columns  1-5  and  read  by 
the  format  F0RMAT(I5).  From  this  point  on,  the  card  order  is 
dependent  upon  whether  or  not  optimization  is  performed. 

f_.  If  optimization  of  model  parameters  is  required,  then  specify 
the  number  of  rainfall-runoff  events,  designated  by  MOBS,  to  be 
used  in  optimization.  The  format  for  reading  it  is  F0RMAT(I5). 

g.  If  optimization  of  model  parameters  is  not  required,  then  model 
parameters  must  be  specified.  Provide  the  lag  parameter  that  is 
denoted  by  PAR.  The  format  for  it  is  F0RMAT(F10.4) . 

h.  Read  the  number  of  rainfall-runoff  events  for  prediction.  This 
is  denoted  by  NOBS.  The  format  for  reading  it  is  F0RMAT(I5). 
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Subroutine  EXOP 


5.  The  purpose  of  this  subroutine  is  to  set  the  stage  if  optimization 
of  model  parameters  is  needed.  The  input  for  this  subroutine  is  given  as 
follows : 

a.  Specify  the  number  of  parameters  for  optimization  denoted  by  N, 
number  of  stage  searches  desired  by  optimization  algorithm 
denoted  by  MST,  control  value  for  printing  of  results  of 
optimization  algorithm  denoted  by  1PT,  convergence  tolerance 
based  on  change  of  objective  function  denoted  by  EPS,  and 
weighting  factor  denoted  by  WF  to  be  used  in  defining  the  objec¬ 
tive  function.  When  IPT  =  0  ,  only  the  final  parameter  values 
are  printed.  When  IPT  =  1  ,  parameter  values  at  each  stage 
search  are  printed.  When  IPT  =  2  ,  parameter  values  at  each 
cycle  search  are  printed.  These  are  read  as  READ(5, .)N,MST, 

IPT, EPS, WF  using  the  format  FORMAT (315 ,F15 . 6 ,F10. 4) . 

b.  Specify  initial  guesses  of  the  parameters  denoted  by  PAR(I), 

1=1, 2,..., n,  where  n  is  the  number  of  parameters  to  be  opti¬ 
mized.  These  are  necessary  to  start  the  optimization  algorithm. 
These  are  read  as  READ(5 , . ) (PAR(I) , 1=1 ,N)  with  the  format 
FORMAT (8F1 0.4) . 

c.  Specify  lower  limits  of  the  parameter  values  denoted  by  PL(I) , 
1=1 , 2 , .  . . ,n .  These  are  ready  as  READ(5, .) (PL(I) ,  1=1, N)  with 
the  format  FORMAT (8F 10. 4) . 

d.  Specify  upper  limits  of  the  parameter  values  denoted  by  PU(I), 
1=1, 2,..., n.  These  are  read  as  READ(5 , . ) (PU(I) , 1=1 ,N)  with  the 
format  FORMAT(8F10.4) . 

6.  The  lower  ana  upper  limits  define  the  range  from  which  optimal 
parameter  values  must  be  derived. 


Subroutine  OBJECT 


7.  This  subroutine  computes  the  objective  function  for  optimization. 
No  input  is  read  in  this  subroutine. 


Subroutine  PRECIP 


8.  This  subroutine  reads  rainfall-runoff  data  for  a  given  watershed. 
Employing  the  information  furnished  by  the  subroutines  NEWTON  and  INFIL,  it 
computes  the  effective  rainfall  and  arranges  it  in  a  proper  manner.  The  input 
to  this  subroutine  is  given  as  follows: 
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a.  Specify  the  number  of  rainfall  readings  in  a  given  event.  This 
is  denoted  by  NNQ.  This  is  read  by  READ(5,.)NNQ  with  the 
format  FORMAT (15). 

b.  Specify  the  date  and  the  watershed  on  whicu  the  rainfall  event 
occurred.  This  is  denoted  by  INF.  The  read  statement  for  this 
is  READ (5, .) (INF (I) ,1=1 ,20)  with  the  format  FORMAT (20A4) . 

c.  Specify  the  volumes  of  rainfall  and  direct  runoff.  These  are 
denoted  respectively  by  RVOL  and  QVOL.  The  read  statement  for 
this  is  READ (5,.) RVOL, QVOL  with  the  format  F0RMAT(2F10. 4) . 

d.  Specify  the  rainfall  hyetograph  where  time  is  given  in  hours 
and  minutes  and  intensity  in  centimetres  per  hour.  Depending 
upon  the  number  of  readings,  this  may  be  specified  on  several 
cards.  The  readings  in  hours,  minutes,  and  intensity  are 
denoted  by  IT1,  IT2,  and  QI,  respectively.  The  read  statement 
for  this  is  READ (5, .) (IT1(I) ,IT2(I) ,QI(I) ,1=1, NNQ)  with  the 
format  FORMAT(4(2I5,F10.4)) . 

e.  Specify  the  number  of  runoff  readings.  This  is  denoted  by  NQQ. 
The  read  statement  for  this  is  READ(5,.)NQQ  with  the  format 
FORMAT (15) . 

_f .  Specify  the  date  and  basin  on  which  the  runoff  event  occurred. 
This  is  denoted  by  INFQ(I).  The  read  statement  for  this  is 
READ(5, .) (INFQ(I) ,1=1,20)  with  the  format  FORMAT ( 20A4) . 

£.  Specify  the  runoff  hydrograph  where  time  is  given  in  hours  and 
minutes  and  discharge  in  centimetres  per  hour.  These  are 
respectively  denoted  by  JTQ1 ,  JTQ2,  and  QOB.  Depending  upon 
the  value  of  NQQ,  these  may  occupy  several  cards.  The  read 
statement  here  is  READ(5, .) (JTQ1(I) , JTQ2(I) ,Q0B(I) ,1=1, NQQ) 
with  the  format  F0RMAT(4(2I5,F10.4)) . 


Subroutine  NEWTON 


9.  The  purpose  of  this  subroutine  is  to  determine  the  Philip  infiltra¬ 
tion  parameter,  sorptivity  S  .  No  input  data  are  specified  in  this 
subroutine . 


Subroutine  BROSEN 


10.  This  subroutine  optimizes  the  parameter  values  for  a  given  set  of 
rainfall-runoff  events.  No  input  is  read  in  this  subroutine. 


Subroutine  INFIL 


11.  This  subroutine  computes  infiltration  using  the  Philip  infiltration 
model.  No  input  data  are  specified  in  this  subroutine. 

Subroutine  XDATA 


12.  The  purpose  of  this  subroutine  is  to  arrange  effective  rainfall 
data  at  equal  time  intervals.  No  input  data  are  required  in  this  subroutine. 

Subroutine  BASIN 


13.  This  subroutine  specifies  and  computes  pertinent  geomorphic  param¬ 
eters.  The  input  in  this  subroutine  is  given  as  follows: 

a.  Specify  the  purpose  of  this  subroutine.  This  is  denoted  by 
PURP  and  occupies  two  cards.  The  read  statement  is 
READ(5, .) (PURP(I) ,1=1,40)  with  the  format  FORMAT (20A4) . 

b.  Specify  general  information  about  the  watershed,  its  location, 
its  type,  etc.  This  is  denoted  by  INF  and  will  occupy  one 
card.  The  read  statement  is  READ(5 , . ) (INF(I) , 1=1 , 20)  with 
FORMAT (20A4) . 

c.  Specify  the  area  and  order  of  the  watershed,  respectively 
denoted  by  A  and  W.  These  are  given  on  one  card.  The  read 
statement  is  READ(5,.)A,W  with  FORMAT(F10.4,I5) . 

d.  Each  channel  element  within  a  watershed  is  identified  by  a 
label  indicating  the  channel  order  and  sequence  number  of  the 
channel  element.  For  example,  1.3  denotes  the  third  channel 
element  of  the  first-order  channel  for  watershed  G  as  shown  in 
Figure  13.  This  identification  of  channel  elements  is  con¬ 
venient  but  not  essential.  Obtain  the  channel  order  having  the 
highest  number  of  channel  elements.  Specify  this  number  of 
elements  by  MAX  and  its  order  of  the  channel  by  OCM  on  the  same 
card.  The  read  statement  is  READ(5 , . )MAX,0CM  with  the  format 
FORMAT(2I5) . 

e.  Specify  the  channel  order  and  the  associated  number  of  ele¬ 
ments,  denoted  respectively  by  OC  and  NC.  Depending  upon  the 
value  of  W,  these  may  occupy  several  cards.  The  read  statement 
is  READ (5 , .) (0C(I) ,NC(I) , 1=1 , W)  with  FORMAT( 1615) . 

fi .  Specify  the  number  of  paths  available  in  the  watershed,  denoted 
by  MS.  The  read  statement  is  READ(5,.)MS  with  F0RMAT(I5). 

g.  Specify  the  path  and  the  number  of  mergers  of  channels  occur¬ 
ring  in  this  path.  These  are  denoted  by  PAT  and  MC.  The  read 
statement  is  READ ( 5 , . ) PAT ( I ) , MC  with  F0RMAT( 1615) . 
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h.  Specify  for  each  path  the  number  of  channels  of  order  i  merging 
into  channels  of  order  j  according  to  the  path  structure.  This 
is  done  by  specifying  Cl,  CJ,  and  ICJ  where  Cl  denotes  the  num¬ 
ber  of  channels  of  order  I  that  will  merge  into  a  channel  of 
order  J  higher  than  I,  and  ICJ  number  of  channels  of  order  I 
merging  into  channels  of  order  J.  Depending  upon  the  number  of 
possible  paths,  this  may  occupy  several  cards.  The  read  state¬ 
ment  for  this  is  READ (5 , . ) (Cl ( J) ,CJ( J) , ICJ (I , J) , J=1 ,MC)  with 
FORMAT (16 15) . 

i.  On  a  card  specify  the  channel  order,  denoted  by  OC.  The  read 
statement  is  READ(5 , . )0C (I)  with  FORMAT ( 15 ) . 

j_.  Specify  the  length  of  each  element  in  a  channel  of  each  order. 
This  is  given  by  NE  and  CL  where  NE  is  the  channel  element  num¬ 
ber  and  CL  the  element  length.  Depending  upon  the  number  of 
channel  elements  and  the  watershed  order,  this  specification 
may  require  several  cards.  The  read  statement  is  READ(5,.) 
(NE(J) ,CL(I,J) ,J=1  ,NCC)  with  the  format  FORMAT (5 (I5.F10.2)) . 

NCC  signifies  the  number  of  channel  elements  of  a  given  order. 

k.  Specify  channel  order,  denoted  by  OC,  on  a  card.  The  read 
statement  is  READ (5 , . )OC(I)  with  FORMAT (15) . 

l.  Specify  channel  element  number  (NE)  and  area  draining  directly 
into  the  channel  (AC).  Depending  on  the  watershed  order  and 
the  number  of  elements,  it  may  take  several  cards  to  make  this 
specification.  The  read  statement  is  READ(5 , . ) (NE(J) ,AC(I , J) , 
J=1 ,NCC)  with  FORMAT (5 (I5.F10.4)) . 

m.  Specify  the  path  number  denoted  by  PAT.  The  read  statement  is 
Read  (5,.)PAT(I)  with  FORMAT (15) . 

n.  Specify  the  path  matrix.  The  spatial  evolution  of  a  water 
particle  through  a  geomorphic  network  of  overland  regions  and 
channels  is  perhaps  best  accounted  for  by  considering  the  over¬ 
land-channel  flow  paths  that  a  water  particle  may  take  from  the 
point  of  its  landing  to  its  arrival  at  the  basin  outlet.  The 
specification  of  these  paths  for  a  watershed  can  be  made  by 
following  the  transition  rules  discussed  previously.  To 
illustrate,  the  overland-channel  flow  paths  for  watershed  G  can 
be  specified  as 
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14.  Here  c<.  is  the  trapping  state.  It  should  be  noted  that  a  water 
particle  always  originates  in  one  of  the  overland  regions.  Furthermore,  a 
water  particle  travels  first  to  the  channel  element  associated  with  that 
overland  region  and  then  continues  its  journey  to  the  outlet  through  higher 
order  channel  elements.  The  last  state  represents  the  trapping  state  as 
exemplified  by  c^  for  watershed  G. 

15.  The  information  on  the  configuration  of  various  overland-channel 

flow  paths  is  supplied  to  the  program  in  the  following  manner.  An  array  con¬ 
sisting  of  r^,  r^.....  r^;  c^,  is  considered.  For  example,  in  case 

of  watershed  G  such  an  array  can  be  written  as 

V  V  r3’  r4;  cl’  C2’  C3*  c4 

A  value  of  1.00  or  0.0  is  inserted  in  place  of  r^  or  c^,  i  =  1,  2,...,  W, 
depending  upon  whether  or  not  r^  or  c^  is  present  in  a  given  path.  If  the 
first  overland-channel  flow  path  for  watershed  G  is  considered,  then  the 
information  pertaining  to  this  path  can  be  coded  as  follows: 

1.0,  0.0,  0.0,  0.0,  1.0,  1.0,  1.0,  1.0 

Likewise,  the  entire  structure  of  overland-channel  flow  paths  can  be  coded  as 


1.0, 

0.0, 

0.0, 

0.0, 

1.0, 

1.0, 

1.0, 

1.0 

1.0, 

0.0, 

0.0, 

0.0, 

1.0, 

0.0, 

1.0, 

1.0 

1.0, 

0.0, 

0.0, 

0.0, 

1.0, 

0.0, 

0.0, 

1.0 

1.0, 

0.0, 

0.0, 

0.0, 

1.0, 

1.0, 

0.0, 

1.0 

0.0, 

1.0, 

0.0, 

0.0, 

0.0, 

1.0, 

1.0, 

1.0 

o 

o 

1.0, 

0.0, 

o 

o 

0.0, 

1.0, 

0.0, 

1.0 

0.0, 

0.0, 

1.0, 

0.0, 

0.0, 

0.0, 

1.0, 

1.0 

o 

O 

0.0, 

0.0, 

1.0, 

o 

o 

o 

o 

0.0, 

1.0 

This  coded  information  on  overland-channel  flow  paths  becomes  input  to  the 
program  mid  c^e^ified  by  I  denoting  the  path  number  and  PATH  denoting  an  array 
corresponding  to  overland  regions  and  channels  appearing  in  the  path.  This 
read  statement  is  READ (5 , . ) (PATH (I , J) , J=1 ,WW)  with  FORMAT ( 1 5F5 . 1 ) . 

Subroutine  LAG 


16.  This  subroutine  computes  the  lag  time  using  a  lag-area  relation¬ 
ship.  No  input  is  read  in  this  subroutine. 
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Subroutine  HOLD 


17.  This  subroutine  computes  holding  times  for  the  paths  available  in 
the  watershed.  No  input  is  read  in  this  subroutine. 

Subroutine  1UH 


18.  This  subroutine  computes  the  instantaneous  unit  hydrograph  (IUH) 
using  the  geomorphologic  formulation.  No  input  is  read  in  this  subroutine. 

Subroutine  CONVOL 

19.  This  subroutine  performs  convolution  of  the  rainfall  excess  with 
the  IUH  to  determine  the  direct  runoff  hydrograph.  No  input  is  read  in  this 
subroutine . 
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APPENDIX  C:  GMHS 


RELEASE  2.0  MAIN  DATE  =  P  3 1  32  23/  1  3/29 

*  **  **********  ************  ****** 

***  MAIN  PROGRAM  GMHS  *** 

******* **************** **«*•••* 

*»*  THIS  PROGRAM  SYNTH FSTZrS  THE  SURFACE  RUNOFF  HYOROGPAPH  USING  *** 
***  THE  MODEL  GMHS  .  *** 

**********  **********  ******************************************** »**•*•** 

01  MENS  ION  DIJPP  t  40)  .0BYI20I  ,08YT(  201  ,  NY  I  2  0)  ,PRt  10) 

COMMON /OBJ /ORY,ORYT,NX 
COMMON / POP /A.MS, W, WW 
COMMON /POM/OT, EX.WF, KCPT 
COMMON /ROSE/PR 
PEAt  KB 
INTFGFP  V,WW 

******************************************************** 

***  FAC  CONVERTS  INCHES  to  CENTIMETERS  *** 

***  CF  CONVERTS  METERS  TO  FEET  *** 

***  PUP  P ( 7  I  SPECIFIES  THE  PURPOSE  OF  THE  PROGRAM  *** 

******  **************************************** ********** 

CE=  3.2808 
FAC  =  2.54 

READ! 5  t  5 )  IPUPPI  I| ,1 =1,40) 

5  FORMAT (20A4I 

WP  I  T E( 6. 10)  (PUPPI IJ ,1=1 ,40) 

10  FORMAT <5X,20A4) 

************************ ********  ************************* 

***  NK  SPECIFIES  THE  NUMBER  OF  BASINS  UNDER  STUDY  *** 

***  DT  SPECIFIES  THE  TIME  INTEPVAL  OF  COMPUTATION  *** 

****************** ************** ****************** ******* 

READI5.35)  DT,NW 
35  FORMAT  (F10. 4,  15) 

WPITEI/,, 45)  dt  ,  NW 

45  FORMAT ( 5X, • TIME  INTERVAL  EOP  COMPUTATION  IN  HOURS  IS  =*,F10.4/5X, 
l'THE  NUMPEP  OE  WATERSHEDS  FOP  STUDY  IS=',T5/) 
on  50  J=1,NW 

A************************* ******  **********  ********************** **  * 

***•  NOP  S  IS  THE  NUMBER  CE  R  A  I N  r  A  LL -P  UNOE  E  FV  ENTS  FOP  A  RASIN  **  * 
******************************************************************* 

PEA n ( 5 . 25  I  NOBS 
25  FOPMAT  j  [5) 

WP!  TO  6,  10  0)  NO  P  S 

100  FORMAT ( 5X,  'NUMBER  Or  R A  I NF AL L-PU NOfF  FVENTS  AVA1LAB1E  ON  THE 
1WAT  EPS  HED  =',  15/) 

******  ****************************  ************************ ******* 

•i**  AA  IS  THE  PHILIP  I  NF  I  L  TP  AT  1 6  N  P  4R  AME  TT  P  APPROXIMATELY  *** 

*•<*  EOUAL  TO  THE  SATURATED  HYDRAULIC  CONDUCTIVITY  ** * 

***  EX=EXPONENT  IN  THE  LAG-AREA  RELATION  *** 

***  NXM  IS  THE  NUMBER  OF  COMPUTATION  TTMF  INTERVALS  **• 

******  ****** ****************** ••**•**•**•* ************ *********** 

R  E  A  D I  5  ,  B  0 )  A  A  ,  E  X  ,  NX  M 
80  FORMAT  I 2F1 0.4, I  5) 

WRITE(6,110)  EX 

110  FORMAT ( 5X, 'EXPONENT  IN  THE  LAG-AREA  RE L AT I 0N= • , F 1 0 .4 / ) 

AA  =  AA*  TAC 
WRITE! 6,  10  5 ) NXM 

105  FORMAT ( 5X, 'NUMPEP  OF  COMPUTATION  STEPS  IS  =',15/1 
WP  ITFt  6, 85 )AA 

85  FOPMATI5X,  'PHILIP  PAPANETEP  A ( CM /HR) =' , F 1 0 .4 / ) 


CL 


onnor>on  nonnoo 


PFLEASF  2.0 


M.\  IN 


DA  T  F 


“1132 


23/13/2° 


RCA0I5  ,  1MKOPT 
15  FflP  MAT  (15) 

CALL  FAS  IN(A,MS,W,  WW  ) 

IF(KOPT.FQ.O)  GO  TO  20 
WPITFI  6  f  65 ) KOP  T 

65  FOPMAT (5X, 'OPTIMIZATION  OF  P  AP  AM  FT  Fn  S  IS  PrQUIPcD  AND  OPTIMIZATION 
ICOOt  IS’, 15) 


***  DEFINE  THF  V  DM  PEP  Or  R AI NF ALL-RUNOFf  EVENTS  TO  0E  USFO  IN 
***  OPTIMIZATION 

***  MCflS  =  N(IMBEP  Or  PAINFALL-°UNOFF  EVENTS  USED  IN  RAP  A  MF  TE  p 
***  OPTIMIZATION 

MtltttttOlIttttttttMtttttttUtXtttttMttttmttMtMMMttM 

REA0(5t25)  MORS 
WRITE! 6*30)  MOBS 

30  FOPMAT ( 5X,  'NUMBER  OF  P A I NF AL L- PU NOFF  EVENTS  FOR  OPTIMIZATION  IS* 


*  *  » 
*** 
*•* 


1.2X, 15/ ) 

CALL  EXOPI MOBS  ,NXM,AA«  A) 


NOR  S=N' OPS— MORS 
GO  TO  70 

20  WP I TE< 6, 75 1KOPT 

75  FORMAT ( 5X, • PAP AMFTERS  APE  KNOWN  AND  NO  nPTIMI?ATTON  IS  NECESSARY'/ 

1  5X  «  ' NO  OPTIMIZATION  CODE  TS',13/1 

-  *  t  *i 

***  PAR  IS  THE  proportionality  FACTOR  AND  FX  THE  EXPONENT  IN  *** 

***  LAG-AREA  P  FL  AT  I  ON  WHFPF  AP  EA  IS  IN  SQUARE  KILO“ETEPS  AND  *** 

***  KOPT  IS  ALSO  USED  AS  A  CONTROL  FOR  OPINTING  O'JTPUT  IN  A  *»* 

***  LAG  IN  HOURS  *** 

***  GIVEN  SUBROUTINE  *** 


READ! 5,55)  PAR 
55  F0PMAT(FI0.4) 

WRITE! 6 ,60) PAP, EX 

60  FORMAT ( 5X»  'PARAMETER  ALPHA  IN  LA  G( HP  5 )  AREA! <"**;>(  P  EL  AT  I  ON  1  S= • 
1  F 1  0 . 4/  5V  , '  F X°ONRN  T  IN  LAG  APEA  RELATION  1S=',F10.4/) 

P  F A  0 ( c , 25JNOPS 
70  CONTINUE 
KOPT=0 


WP  I  TH  6,  96  )MORS 

96  FORMAT ( 5X,  'NUMBER  OF  PAINE  ALL  RUNOFF  EVFNTS 
115/  ) 

PAP  =  PP  f  I  ) 


FOR  OPFOIFTICN  IS'. 


CALL  L  AM  A  ,PAR  ,EX,  KP  ) 

CALL  HOLD! MS, W, A, WW, KB, KOPT) 

CALL  I  UHl M S , WW ,OT , KOPT  j 
00  40  1=1, NOBS 

CALL  PPEC I  PI NNQ,DT,AA,OBQ,OBOT fA ) 

CALL  X  DAT  A  (  NNQ  »NXM  ,f)T  I 

CALL  CONVOL(DT,NXM,QP,QPTf KOPT.A ) 

0FPp=i obo-op i/nno 

QTERP=(OPQT-QPT)/OBQT 
WP ITr( 6,115)  OBQ , QP, QEPP 

115  FORMAT  (  5X,  'OfVS.PFAK  01  SCH  A  °G  F  I  CM  /HR  )  =  •  ,  T  1  ">  .  4 , 2  X  ,  •  C  OMP.  P^AK 
IDISCHAPGE= • .F 10.4, 2X, 'RELATIVE  F  PF OR  =  '  , F 1 o  .4 / ) 

WPITFI 6,  120)  ORQT, QP  T , OTFR  P 

120  FOPMATISX,  'OP.S.  PEAK  T  I  MR  I  M  t  N )  =•  ,F  1 0  .  3 , 2  X  ,  •  COMP.  PFAK  T  I  MR  (MINI  = 
1 *,F 1 0. 3, 2X, 'PEL  ATI VE  ERROR  =',F10.4Z) 

40  CONTINUE 


RELEASE  2.0  MAIN  r^TF  =  8313? 

50  CONTINUE 
STOP 
END 


C.3 


22/47/28 


nnnnn  riooortr'rtnrinri  nnnn 


PELEASF  2.0 


exop 


DATE 


8X132 


23/  13/?*} 


SUBROUTINE  EXnPC'OBS.NXM.AA,  A) 

******************  ******************************  ****** 

***  THE  PUPPnSE  OF  THIS  SUBROUTINE  IS  TO  PROVIDE  PFPTJNFNT  *** 
***  INFORMATION  REQUIRED  BY  TH  F  OPT  I  Ml  7  AT  I  "N  ALGORITHM.  *** 

********************  ********** ******************************  ****** 

DI MENS  ION  Til  100) , QI  (  1  00  5  , XI  2000  1, PI  10,10001  ,PAP« 101 , PL  I  10 1 , PUI 1 01 
1,  NX  I  20 1 ,OBY! 201.0BYTI 201 
COMMON/VS/TI.OI.P 
COMMO;  /PPM /DT , EX.KF, KOPT 
COMMON / OBJ /O BY ,OBYT , NX 
COM  MON /PAR  A/PAP, PL ,pu 
COMMON /VSX/X 
INTEGER  W,  UN¬ 
READ!  5,5  IN,  M$T,  I  PT,E  PS.  WF 
5  FORMAT ( 3I5,F15.6»F 10 .41 

********************************  **  **************  ****  **************  ** 


***  N=NUMBER  OF  PARAMETERS  *** 

***  MST=  NUMBER  C,F  STAGE  SFARCHES  DESIRFD  *** 

***  1PT=  0 - ONLY  THE  FTNAI  PARAMETER  VALUFS  PRINTED  *** 

***  I PT=  1 - INTERMEDIATE  VALUFS  OF  EACH  STAGF  SEARCH  PRINTED  *** 

**»  IPT= 2 - INTERMEDIATE  VALUES  OF  EACH  C YC L F  SFAPCH  PRINTED  *»* 

***  EPS= CONVERGENCE  TOLERANCE  BASED  ON  THE  CHANGE  OF  THE  *** 

***  OBJECTIVE  FUNCTION  *** 

***  WF=  WEIGHTING  FACTOR  TO  BE  USED  IN  DEFINING  OBJECTIVE  •** 

***  FUNCTION  *** 


****************************** ******  ********************** ********** 

WRITFI 6,10  IN, Fps 

10  FORMATC5X, 'NUMBFP  OF  PARAMETERS  =• , l 5, 3X , • CCNVERGENC E  TOLERANCE 
IllMIT= •, El  5.7/1 
WRITE! 6,15)MST,IPT 

15  FORMAT (5X, 'NUMBER  OF  STAGF  SEAPCHES  S^EC IF IED= • , I5.5X, • CODE  FOR 
IPR l NT I NG  PAPAMETEP  VAL  UF  S=  •  ,  1 5 /) 

WRITC( 6,45  )KF 

45  FORMAT (5X, 'WFIGHTING  FACTOR  USED  IN  DEFINITION  PE  OBJECTIVE 
1 FUNCT I ON= ' , FID. 4/1 

************************************************ 

***  PAP( I  )  =  INITIAl  GUFSS  OF  I-TH  PARAMETER  *** 

***  P  L ( I  1 =  L  OW  F  R  BOUND  OF  I-TH  P  AR  A" FT F P  *** 

***  P U(I1=UPPEP  BOUND  OF  I-TH  PAPAMETPP  *** 

******************  ****************************** 


20 


25 

30 

35 


37 


PFADI5 ,20) (PAP  I  I),  1= 1,N) 

FORMAT (PFIO. 4) 

READI5 , 20)  (PLm,!=l,N) 

P  EADI  5 , 20 1  IPUin,I*l,N) 

WPJTFI 6,25)  (PARI  I  1, 1=1, N) 

FORMAT  I 5X , •  INITIAL  PARAMETER  GUESSFS  ARE • ,P F 1 0 . 4/ 1 
WRITEI6.30)  I  PL  I  1 1 , 1 =1 ,N 1 

FOP  MAT  I 5X ,  •  LOWER  LIMITS  OF  PARAMETER  VALUES  ARE • , 8F LO .4/ 1 
WRITFI  6,35)  (PUin.I=l,M) 

FORMAT (5X, 'UPPER  LIMITS  OF  PARAMETER  VALUES  ARE',  8F10.4/1 
DO  40  1=1. MOPS 

CALL  PRECl P(NNO,OT,AA,OBO,OBQY,A) 

OBY  I  n  =  OBO 

OBYT I  1  )  =  OB0T 

CALL  X  DA  TA l NN0 , NXM  »DT 1 

NX  (  I  1  =  NX M 

on  37  J  =  1 , NXM 

PI  I , J)=X(J  1 


C4 


REL  EASE  2.0 


CXCIP 


DATE 


8313? 


22/47/28 


40  CONTINUE 

CALL  BROSENIN.MST, I P T, EP S, MOOS  I 

RETUPN 

END 


C5 


nononnn 


RELEASE  2.3 


OBJECT 


DA  TC 


R31  3? 


23/ 13/29 


20 


19 


30 

25 


35 


SURPOUT INC  OBJECT! VALUE, MORS,  N I 

************<<<.****  **********************************************  ******** 

***  THIS  PPPGPAM  SPECIFIES  THE  OBJECTIVE  FUNCTION  TO  BE  USFO  IN  *** 

***  OPTIMIZATION  OF  PARAMETERS  *** 

***  THF  OBJECTIVE  FUNCTION  IS  DEFINED  HERE  AS  THE  SUM  OF  SOUARFS  *** 

***  OF  DEVIATIONS  BETWEEN  OBSERVED  AND  C0MPARED  OFAKS  AND/OP  *** 

**«  THEIP  TIMES  OF  OCCUPENCF  *** 

************************  **********************  ******** 
01  Mr MS  ION  op  I  10) , NX  I  2  0) ,P{  1 0,  10  00) ,PPEDY(20 ) ,X  t  2000 l,ORY< 20) 

1 .OBYTI  20).  PPFDTI 20), TI { 1001 ,01 (100) 

COMMOM/ROSE/PR 
COMMON /OBJ /OBY .OBYT.NX 
C OM  MON  /POP  /  A  ,  M  S  ,  K,  WW 
COMMON /POM /OT, EX.WF, KOPT 
COMMON /VS/T1 , Q I , P 
COMMON /VSX/X 
REAL  KP.KB 
INTEGEP  W.WW 

WPITEI  6,5)  (PPII), 1=1, M) 

F0PMATI5X,  'INITIAL  PAPAMETER  GUESSES  APE  •  ,8F  10 .4/ ) 

KP=  PPI  1) 

CALL  L AG( A.KP, FX.KBI 

CALL  HOLOl  MS.W, A.WW.KR .KOPT) 

CALL  IUHIMS.WW.OT.KOPT) 

00  15  1=1, MOBS 
NXM  =  NX ( I  I 
DO  20  J= 1 , NXM 
X  <  J  1  =P  *  I  .  J  ) 

CALL  C  ONVOL (nT .nxm.op.qpt.kopt ,A ) 

PREOYI I )=QP 
PREDT (  I  >=QPT 
CONTINUE 
WP I TF( 6, 10) 

FOPuAT  t  5X»  'OBSERVED  PEAK* , ?X , • PP ED K TFD  r>f AK  •,? X, • FP ROR= CO-OP' ,2X 
1,'CRS.PFAK  TIME*  ,2X, 'PRED.  PEAK  T1M*'»2X»*EP  ROR=  '  PO— T°C.  •  / ) 

VALUE  1 =0.0  0 
VALUE?=0.03 
DP  25  1=1, MO^S 
DEVO=CRY<  I  I-oc-EDYI  I  ) 

TEMP=DEVP*DEVO 

vai  uri  =valuei+tfmi> 

OEVT=OBYT(  I  I  - D P. E 0 T (  I  I 

tempt=dfvt*devt 

VALUE2=VALUE2+TEMPT 

WPITEI 6, 30  I  OBY I  I ),PPEDY( I ) ,OEVC, OB YT (I) , »P EDTl I ) ,DEVT 
FORMAT I6I5X.F10. 4) ) 

CONTINUE 

VALUF=WF*VALUF1*(1 . O-WF ) *VAL UE 2 
VMQ=VA  LUE1 /FLOAT  I  MOB  S) 

VMT=VALUE2/FLnATIM0RSl 
WRITE! 6, 35) VAI UF.VMO.VMT 

FOP  MAT ( 5X,  '  OBJECTIVE  FUNC TI ON  V AL UF  =  ' ,E 1 5 .7/5 X , •  MEAN  SQUARE 
1ERR0P  IM  HYDPOGR APH  PFAK  =  * ,  F15.7/5X,'  MEAN  SOUARF  ERROR  IN  TIME 
2T0  HYDPOGPAPH  PFAK=',  E 1 5. 7/  ) 

RETUPN 

END 


C6 


nnn  non  noon  oonooon  noon  noon 


RELEASE  2.0 


BASIN 


DATE 


R  3 1  3? 


23/13/29 


SUBROU 

****** 

*** 

*** 

****** 

DIMFNS 

I  AC ( 20, 
1 ,  PA  T  (2 

COMMON 
COMMON 
1NTFGE 
COEF  =  1 
CFF= 1 . 
****** 
*** 

*** 

****** 

RFA  0( 5 
5  FORMAT 
WPITEI 
READI5 
10  FOP MAT 
HP  IT  FI 
15  FORMAT 
****** 
*** 

*** 

*** 

*** 

*** 

****** 

REA0I5 

20  FORMAT 
A=A*CO 
WP  IT  FI 

21  FORMAT 
****** 
*** 

*** 
****** 
R  F  AD ( 5 

22  FORMAT 
WP  I  T  FI 

25  FORMAT 

II  5/ I 
****** 
*** 

****** 

RE  ADI  5 
30  FORMAT 
WPITEI 
35  FORMAT 

****** 

*** 

****** 


TINE  MAS IMA, MS, W,WW> 

********** ********************************** ************** 

THE  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  POOVtDF  PFPTINFNT 
GEC-MORPHOLOG IC  INFORMATION  PFQUIRFD  BY  THE  GMHS  MO^FE 
**************  ******************************** ************ 
I ON  PUP P (401 , I NF120I ,NC ( 20  I, fC  J<  20,2  0)  ,CL( 20,50) , ALI 20  )  , 
501,  SAI 201, ARI  20) ,PC (2  0,20 ) ,  MP (201 ,PSI 20) , PA (2  0) , PAT HI  20,2 
0) ,F II 20), CJ (20) ,OC (20) ,NE (23) 

/ONE  /AR  ,  Nr.  ,  AL 
/TVD/PATH, PS 
R  W,WW,CI,CJ,PAT,OCM,OC. 

.0 

0 

****************************************************** 
PUPP(I)  SPECIFIES  THE  PUPPOSE  OF  THIS  SUBPOUT  I NF  *** 

INF  I  I )  GIVES  NAME  AMD  LOCATION  OF  THE  "ASIN  **» 

********************** ************ ******************** 

,5)  (PUPPIII ,1*1,431 
I 20A4I 

6,5)  (PllPPm.  1*1,40) 

,10)  (IN FI  I), 1*1, 20) 
l  2  3A  4  I 

6,15)1  IMF! 11,1=1,20) 

I 5X , 20A4  ) 

************************  **  ******************************** 

THE  LENGTH  AND  AREA  APE  IN  KILOMETERS  AND  SQUARF  K  I L  0  M  FT  F° 
RESPECTIVELY.! F  NOT  THEN  CHANGE  CFF  AND  COEF  FROM  1.0  TO 
APPROXIMATE  VALUES 

A  DENOTES  BASIN  AREA.  IE  A  RFA  IS  NOT  IN  SQUAPE  K I LOM  FT  E°  S 
THEN  CONVERT  IT  TO  THESE  UNITS 

************************  **  **  ****************************** 

,20)  A , w 

I F 10.4,  I  5) 

E  F 

6,21)  A , W 

(5X,  •BASIN  A  RE  A I  SO  K M I = •  ,F  1  0 . 3 , 5 X , • CRD EP  OF  THE  BASIN=',I5 
*<■  ************************************  **  ******  ************ 

OC  '*  GIVES  THE  CHANNEL  ORDER  WHICH  HAS  H'GHEST  NUMR  E®  OF  xh 
ELEMENTS  DENOTED  BY  MAX 

***£**««**»****•*•*• ******  ******************************** 

,221  MtX.OCM 
(214) 

6,25)  MAX, OC  M 

I5X,  'MAXIMUM  NUMBER  OF  E LE MC NT S= • , 1 5 ,2 X , ' I N  CHANNEL  OPDER= 


******* 
S  *** 
*** 
*** 
**  • 
*** 

******* 


/) 

******* 
r  *** 
*** 

**  ***** 


************ 
NC  IS  NUMBFP 

************ 
,30)  (or. (i), 
( 16151 

6,35)  IOCII) 
I 5X, 'CHANNEL 
************ 
COMPUTE  THE 

************ 


**************  ************************ ********* 

OF  CHANNEL  ELEMENTS  AND  OC  CHANNEL  ORDE P  *** 

****** ****************************** *********** 
NCI  I  ),I  =  1,W) 

,NC( I) ,1  =  1  ,W) 

OROER=', 15, 5X, 'NUMBER  OF  CHANNEL  EL EMENTS  =  ' , I  5) 

**************  ******************** ************** 
NUMB  EP  OF  POSSIBLE  PATHS  IN  THF  WATERSHEDS  *** 

**************  ****** ******************  ******** •• 


MS=2** (W-l ) 

WPITEI  6 ,  5 5 ) M S 

55  FORMAT  I 5X, 'NUMBER  OF  POSSIBLE  PATHS  IS=',I5/) 


C7 


nnnnno  <*>nn  nnnn  nnononmon  r>nn 


PELFASE  2.0 


PA  SIN 


DATC 


83132 


23/ 1 3/29 


M«»**»nt**t»«**»  **************** *•***•*•••***••****•••*****•** 

***  READ  THF  NUMBER  OF  AVAIIAPIF  PATHS  IN  THT  WATERSHEDS  *** 

************** ************************************ ************** 

READ! 5  « 133)  MS 
130  FERMAT (15) 

WPITFf 6, 135)  MS 

135  FnpMAT ( 5X, 'NUMBER  OF  AVAILABLE  PATHS  IN  THIS  WATERSHED  IS  =',15/1 

********************************  ********************************  ******** 

***  SPFCIFY  THF  MUMPER  OF  CHANNELS  OF  OROCR  I  “ERGING  INTO  CHANNF  *** 


***  LS  OF  ORDEO  J  ACCORDING  T3  THF  ACTUAL  °A TH  COMPOSITION  *** 
***  PAT  DENOTES  PATH  NUMRE P.  *** 
***  MC  DENOTES  NUMBER  Or  MFRGEPS  *** 
***  Cl  AND  CJ  DENOTE  ORDER  OF  CHANNELS  I  MERGING  INTO  ORDER  J  *>** 
***  OF  CHANNELS  J.  ICJ  IS  THE  NUMBER  OF  CHANNELS  MFPGING  *** 
***  PC  IS  THE  PROPORTION  OF  CHANNFLS  OF  ORDFR  I  MERGING  INTO  *** 
**♦  CHANNELS  OF  OPDEP  J  *** 


**************************************************************** **** **** 

00  36  1=  l, MS 
READI5  »  3T 1  PAT(I),MC 

RE ADI  5, 37)  ( C II J ) , CJ I J ) , IC J{ I , J)  , J  =  1 ,MC ) 

ME!  I )  =  MC 

37  FORMAT ( 161 5) 

WP I T FI  6,  39)  PATl  I),MC 

39^FPRMAT I5X»'THIS  PATH  IS  =  • , 2X  ,1 5. 2 X , • NUMBER  OF  MERGERS  IS',2X,I5 

HR ITEI  6, 38 )  (C  1 1 J)  ,CJ<  J) , ICJ< I ,J ), J=1,MC) 

******  ************ ******** ********************************** ******** 

***  COMPUTE  PROPORTION  OF  CHANNFLS  OF  ORDFR  I  MERGING  IN  TO  *** 

***  CHANNEL  F  S  OF  ORDER  J  ACCORDING  to  PATH  STRUCTURE  *** 

********************************  **********  ******** ****************** 

DO  95  J  =  l ,  MC. 

JJ  =  C  I(  J  I 

PCII,  )=  FLOAT ( ICJ II  ,J) ) /FLOAT  INC! JJ ) ) 

95  CONTINUE 

WRITFI  6, 10 D) (CI( J) ,CJ( J) .PCI  I, J)  ,J=1 ,MC) 

100  FORMAT ( 5X, 'PROPORTION  OF  CHANNELS  OF  non EP'*2X,!5»2X.' MERGING  INTO 
1  CHANNF1S  OF  OROEP 2X, I 5,2X, • IS=' ,C1C.4/) 

36  CONTINUE 

38  tOPMATI 5X,  'NUMBER  nF  CHANNELS  CT  OPDEP ■ , 1 X ,  I  2 , 2 X , • ME RG I NG  INTO  CHA 
1NNELS  nr  OPOFP' ,1X»I2»IX,' IS='.2X, 15/) 

********************************  ******* 

***  SPECIFY  THF  CHANNEL  LENGTHS  *** 

*************************************** 

00  41  1=1,W 

NCC=NC ( I ) 

REA0I5.880)  OCIII 
880  FORMAT  I  1 5) 

************************************************************************ 


***  NE  IS  THE  SEQUENCE  NUMBER  ASSIGNFO  TO  A  CHANNEL  ELEMENT  OF  A  *** 
*•*  GIVFN  ORDER  AND  CL  THE  LENGTH  OF  THIS  FLEMFNT  *** 
***  IF  CL  IS  NOT  SPECIFIED  IN  KI LOME  TF  RS  THEN  CONVFPT  IT  TO  THESE  *** 
***  UNITS  *** 


************************************  ******  ************************  ****** 

READ (5, 42)  CNF  I Jl.CLI I  ,  J  )  ,  J=  1  ,  NC  C I 

42  FORMAT  I 5( I 5.F10.2)  > 

WP ITFI 6,43 )  OC (I ) 

43  FOP  MAT  (  5X »  '  CHANNFL  OPDEP.  IS  =',2X,I5/) 

DO  140  J=  1 »  NCC 


C8 


uuc  uuuoouu  uuu  cuuuu 


PELEASr  2. n 


BASIN 


OATE 


A3  l  32 


23/ 13/29 


140  CL  (  I,J  l=CL(  I, J  >*CFF 

HR  IT  FI  6,44  )  I 'I  FIJI  »CLI  I  »  J )  »J  =  1  »NC.f  ) 

44  FORMAT ( 5X, 'CHANNEL  ELEMENT  NUMBE R= * , 12 ,2 X, 

1 'ELEMENT  L  ENGTHI KM )= • ,F10. 3/ I 
41  CONTINUE 

****** ****************************** ***9******************** 

***  COMPUTE  AVERAGE  LENGTH  OF  CHANNELS  OF  EACH  ORDER  *** 

no  46  1=1, w 
SUM=0.0 
NCC=NC (  I  ) 
on  49  J=  1,  NCC 
SUM  =  SU  M+CL I  1 , J  ) 

49  CPMTINUF 

************************************************************************ 


***  Al  IS  THE  AVCRAGC  LENGTH  OF  CHANNELS  OF  A  GIVEN  ORDER  *** 
***  AC  IS  THE  DRAINAGE  AREA  OF  A  CHANNEL  ELEMENT. IF  THE  AREA  IS  *** 
***  NOT  IN  KILOMETERS  THFN  CONVERT  IT  TO  THESE  UNITS  *** 
***  SAU)  IS  THE  OVERLAND  AP  EA  OF  ORDER  I  ♦** 
***  APIII  IS  THE  RATIO  OF  OVERLAND  AREA  OF  ORDER  I  TO  BASIN  AREA  *** 


**************************  ******  **************************************** 

Al  (  I  1=  SlIM/ NCC 

46  CONTINUE 

WRITE! 6,471  (  I , AL!  1 1  , I =1 , W » 

47  FORMAT < 5X,  'CHANNEL  ORDER  = « , 2  * , 1 5 ,2 X, ' A VE R A GE  L ENGTH= • , F 1 0. 3/ I 

******  ************************************  ****** 

***  SPECIFY  AP  EA  OF  EACH  CHANNEL  ELEMENT  ♦** 

******  **************************  **************** 

DO  50  1  =  1,  W 
NCC  =  NC 1  I ) 

P.FADI5, 8R0)  OC  (  I  I 

P FAD (5, 401  i NEIJ I, AC  It ,JI , J=l, NCC) 

40  FORMAT ( 5! I  5, FLO. 41 1 
no  145  j=1.ncc 
145  Ar(  | , J  )  =  AC ( I , J )*COEF 
WRITE! 6, 45 )nc!  1  1 

45  FOR  M  AT ( 5  X ,  • CH ANN  El  CROFR  IS  =«,I5/I 
WRITE! 6, 4P I  (NE(JI,AC!I,J),J=1,NCC) 

68  FORMAT ! 5X, 'CHANNEL  ELEMENT  NUMBER  =',2X,!5,2X, 

I 5X , ' AP  EA  (  SO  KMJ  =  * ,r 10.3/) 

50  CONTINUF 

********************************  ********* 

***  COMPUTE  SURFACE  ARFAS  OF  CHANNELS  OF  EACH  ORDFR  AND  THEREFORE  •** 
***  OF  EACH  PATH  *** 

***  COMPUTE  OVERLAND  AREA  OF  EACH  PATH  *<■* 

************************************************************************* 

DO  56  1  =  1, W 
SUM=0. 0 
NCC  =  NC (  I  ) 

DO  57  J=l, NCC 
SUM=SUM+AC ( I , J  ) 

57  CONTINUE 
S  A!  I  )=  SUM 
API  11=  SAM  I/A 
56  CONTINUE 

WRITE! 6,R5 )  I  SA(I) ,OC(tl ,AR( I) ,1=1 ,W) 

86  rrpMAT ( 5X, » AREA  OE  PEGIONISQ  KM) =» ,F 1 0 . 2 ,?X , » OF  OPDE R= • , 15 , 2X, • P AT 
110  or  REGION  ARFA  TO  BASIN  AREA  =',F10.4/I 


C9 


nnnn  r>r>o 


RELEASE  2.0 


BASIN 


DA  TF 


B  3  1  32 


23/13/29 


*********************************** 

***  SPEC  IFY  THF  PATH  MATRIX  ♦»* 

*********************************** 

WW  =  2*W 
WRITE! 6,65  1 

65  FORMAT ( 5X,  'THE  PATH  MATRIX  IS  AS  F0LL0WS*/1 
OF!  7  5  I  =  1 ,  M  S 

********************************************************************* 

***  PATH  IS  THE  PATH  MATRIX  WHOSE  DIMENSIONS  WOULD  RE  NUMBER  OF  *** 
***  PATHS  X  TWOFOLD  basin  ordfr  *** 

********************************************************************* 
PEADI5,61)  P  AT ( I  1 
61  FOP  MAT (  15) 

WPITTI  6,71  I  PAT<  I) 

71  FOP  MAT (5X,  'PAT H= ' ,  12/) 

PFADI5.6D)  IPATHI I, J> ,J=1  ,WW) 

60  FORMAT  I  15F  5.1) 

WPITEI  6,70  1  (  PATHU  ,J|  ,  J=  l,WW) 

70  FORMAT! 5X, 10F10. 1) 

75  CONTINUE 

DO  120  1=1, MS 
DO  125  J=1,W 

1FIPATHI I, Jl.LE.O.O)  GO  TO  125 
PA  I  I  |=  AP  (  J  ) 

125  CONTINUF 
120  CONTINUE 

00  105  1=1, MS 
TFMP=1 .0 
MC=MFI  I  1 
OP  110  J  =  1  »  MC 
PP0D=PC(  I, J  >»TEMP 
TFMp=pPOD 
110  CONTINUE 

PS(  I )  =  TFMP*PA(  I ) 

WPITFt  6, 11 5)  I.PSI  II 

115  FORMAT ( 5X,  'P  ATH  NUMB FP = • , I  5 , 5X , *  PA TH  PRPB A BI l I T Y* >  , F  10 . 4/ I 

105  CONTINUF 
RETURN 
END 


1 


CIO 


uuu  uo-> 


RELEASE  2.0 


LAG 


DATE 


8313? 


23/ 13/?o 


SUBROUTINE  L AG< A ,P AR  ,EX , KB  ) 

******************************************************* 
***  THIS  SUBROUTINE  COMPUTFS  THF  BASIN  LAG  T I MF  *** 
******  **************  ************  **  ********************* 
REAL  KB 

********************************************* 

***  KB  IS  THF  BASIN  LAG  TIME  IN  HOURS  *** 
********************************************* 
kb*pap*a**ex 

HP  ITFf  6,5) A, KB 

5  FORMAT (5X,  'BASIN  AREA  IN  SQ  KM  = • , FI  0.3 , 2X , • PASI N  LAG  IN 
IF  10 . 2/  1 
RETURN 
END 


HOURS 


Cll 


mn  ,*nnn  noon  rio.^o 


REL  FAS  f  2,0 


HOI  0 


DATE 


P  31  32 


23/13/29 


SUBROUTINE  HOLD!  MS  ,  W ,  A  ,K  W ,  KB  ,KOPT) 

<:***** ****************4************************************************^** 

***  this  suqpoutine  computes  holding  times  of  ovfpland  and  channel  *** 

***  FLOv:  PEG  IONS  **< 

******************  ******************  **********  ********  *********< 


D  I  MENS 
IPS!  20) 
I  NT  EOF 
PFAL  K 
COMMON 
COMMON 
****** 


ION  KR(?0) ,KC( 20), AL (20) , NCI 20),KKI20I . A R( 20 >, PATH (2 0,20) , 

P  K.KW 
c  *KT  ,KK,KP 
/CNF / AP ,NC,AL 
/T^O/P ATH.PS ,KK 

*************************************************************** 


***  SFT  UP  LOGIC  FOP  COMPUTING  HOLDING  TIMES  FOP  C V  FPL  AN  D  AND  *** 
#**  CHANNEL  PEGIONS  *** 

ft *************** ********************************************** ******* 


DO  15  J=l, W 

KPI  J  1=  (  API  J  )«=A/!2.0*NC  (  J)*ALI  J  ))  >**(  1.  0/3.  0) 

KC I  J  PH  )  =  AH  J  )**(  1.0/3.01 
15  CONTINUE 

DO  2P  1  =  1, WW 
IF!  I  .CT.K)  GO  TO  25 
KM  I  l=KP  (I  1 
GO  TO  23 
25  KK(  I )  =  KC (  I  | 

20  CONTINUF 
SSUM=0 . 0 
on  5  1=1, MS 
SU«=0. 0 
DO  10  J=1,WK 
PK=PATH(  I ,  J  )  *KK(  J) 

SUM  =  SUM*-PK 
10  CONTINUF 

PT=SUM»PS(  I  I 
SSUM=S  SUM+PT 
5  CONTINUF 

******************************************** t****** ********** ********* 

***  DETEPMINE  THF  COEFFICIENT  AG  APPEARING  IN  T  H  P  HOLDING  TIME  *** 
***  P FLAT  ION  SHIP  •  »* 

AG  =  K  5/ SSUM 

****************  **  **  ******  4***4*  **  ********  *******  ***************  ****** 

***  DETER M INF  THT  HOLDING  TIMES  OF  OVEPLANO  AND  CHANNEL  PEGIONS  »*« 

****** ****  ********** ****** **** ** ************************************** 

DO  33  1  =  1,  Vi w 
TEMP=l .0/( KK (  I  )  *  AG  1 
KM  I  )  =  TEMP 
30  CONTINUE 

IF(KOPT.NF.D)GO  TO  53 
WPITFI 6,35)1  I , KK< I ), 1  =  1, W) 

35  F0PMATI5X, 'OVERLAND  REGION' , IX, l 5, 5X , 'HOLDING  TIME(HOURS)  IS=',F10 
1.2/) 

DO  60  1=1, W 

KP  I T  E  (  6 , 40  )  I  ,  KK ( I +W  » 

60  CONTINUF 

40  FORMAT (5X,  'CHANNEL  OROER  I S= • , 2X , I  5, 2X , • HOLD  1 NG  TIME  HOURS  IS  =', 

1  F  1 0 .2/  1 
50  CONTINUE 
RFTUPN 


C12 


PELEASF  2.0 
END 


HP  ID 


DATF  =  83132 


22/47/28 


C13 


onn  ooino  nooo 


PELEASC  2 


I  UH 


DATF 


P3132 


2  3/ 13/29 


oUPROUTINF  IUHIMS.WW  ,DT,KOPT) 

***  THIS  SU«PPUTJNE  COMPUTES  THE  INSTANTANEOUS  UNIT  HYDRQGRAPH  *** 

***  USING  GEOMOPPHIC  FORMULATION  *** 

***  *  **  **  **  **  **  ***  $  ********  *  *************  **  *****  *  ******  4c  *  **  **  **  **  **  **  *  * 

niMENS I nr  PATH! 20, 20 >,KK I? 0),NK(  20 >,f (20,20) ,P S I  20 ) , HI  2000 ) 

1 ,PK (  20 ,201  , TT(  1000  ) 

0  F  A 1  KK 

COMMON /TWO /PATH, PS ,KK 
COMMON /THRFF/H.NH 
IMTEGEP  WK 

«««  DFTFPMINE  THE  COEFFICIENTS  C  (  I  ,  J 1  FOP  FACH  PATH.  THIS  RFOUIPFS  *** 
***  AMOTHFP  DETEPMI NATION  OF  M  AT  R I  X  OT  F  Of F I  Cl EMTS  ASSOCIATED  WITH  **< 
***  EACH  PATH. IN  THIS  MATRIX  ONLY  NONZERO  VALUES  VILE  BE  PFTAINFD  * ** 

******  ********  **  ******  ttttti 

KNUMR= 1000 
00  5  1=1, MS 
I  J=  0 

00  10  J  =  l,w*»‘ 

iF(  PATH!  I  ,  J  I  .LE.  0.  0)  GO  TO  10 
I  J=  1  J«-  1 

PK (  1,1  J )=  KK( J  | 

MR (  I  1=  I  J 
10  CONTINUE 

S  CONTINUE 

DO  IS  I  =  1 ,  M S 
NKK  =NK  I  I  ) 

TFMP=1 .3 
DO  PS  J  =  1 , NKK 
PPOD=TEMP*PK( I , J  * 

TfM  P=P  POD 
B5  CONTINUE 

DC  23  IK=l,NKK 
STOP  =  1.0 
DO  2S  J= 1 , NKK 
IK  J.EO.IK  )  GO  TO  25 
0EN  =  STro*I PKt I ,J  )-PK(  I  ,  I  K  I  I 
STOB=PFN 
2S  CONTINUE 

C(  I , IK  I  =  PF  00 / S  TOR 
20  CONTINUE 
15  CONTINUE 

*************************** 

***  COMPUTF  THE  IUH  *** 

**********  ***************** 

U=  l 
HMX=0. 0 
T  =  0. 1 
HI  I  J  )  =  0 . 0 
TT ( 11=0.0 
45  CONT INUE 
SSU  M=D . 0 
00  30  1=  1, MS 
NKK  =  NK  (  I  I 
SUN=0. 0 
DO  t 5  J= 1 , NKK 

PPrD=C (I, JIT  EXPI-PKI 1 , JI*T)*PS  ( II 


uuu 


RELEASE  2.0 


!  UH 


DATE 


8  3132 


23/I3/?o 


sum  =  sum+pp  no 

3s  continuf 

SSUM=S  SUM«-SU‘-1 
30  CONTINUE 
I J=I J+  l 
H(  I  J )=  SSUM 
TTI ! J) =  T 

IE(  HMX .IT.  H<  [  J  I  I  HM  X=H!  I J 1 
HM.H  =  0.  01*HMX 

IF!  HIT  J)  .CE.KNtIMB)  00  TO  40 
!E( HI !  J ) .L  T.HMN I  GO  TO  40 
T=T +0T 
r,n  to  45 
40  CONTINUE 

IF (  IJ.LT. 10)  GO  TO  33 
NH=  I  J 

I F| KOPT.NE.OJGO  TO  00 
00  1  50  I  -  1 ,  MS 
NKK=NK I  I  ) 

WR  I  T  E(  6,  60  )  1 

60  FORMAT  1 5X, 'THIS  PATH  IS  *,15/) 

WRITE(6,80I  <P<(  ! , J ) , J  =  1 ,NKK  ) 

80  F0PMAT(5X,<  PATH  HOLDING  TIMES  A  RF  • /5 X , 1  OF 1 0 . 4/ ) 

150  C  ON  TIN  UE 

DO  155  1=1, MS 
NKK=NK ( I ) 

WRITE!  6,65)  I 

65  FORMAT  I 5X,  'THIS  PATH  IS  ',15/1 
WRITE!  6,70)  ( C  (  I , J  ) , J= 1 , NKK ) 

70  FPPMnT (5X,  'PATH  COEFFICIENTS  CII.J)  ARE  * / 5X , 5 E20. 7/ ) 

155  CONTINUE 

WP I TE(  6, 120) 

120  FORMAT ( 5X, 'THE  I N S TA NTAN FPUS  UNIT  HYDROGPAPH  IS  AS  FOLLOWS*//) 
***  CONVERT  TIME  TO  MINUTES  *** 


On  125  1=1, NH 
125  T  T (  I  )  =  TT( I  ) *  50 . 0 
WP  ITE!  4,  130) 

1  3P  FORM  AT ( 7X,  *T  IME* ,7X,  *1UH* , PX  ,*  TI ME • ,7X, • IUH*  ,RX, 'TIME'  ,7X, * IUH',  7 

IX,  '  T  IME  '  ,8  X.  «I'JH'  ,  7X  ,  *  TIME  •  ,9X,‘  !UH*  /) 

WRITE! 6, 135) 

1  35  FORMAT ( 8X,  »MN», 7X,  *1 /HP  * , 8X, *MN*  ,8X, *1 /HP'  ,8X, *MN* ,»X, • 1 /HR*  ,7X,  • 
1MN '  , 8X ,< I/HP ',  7X,  *MN  •  .OX, ' 1/HR  ■/  i 
WRITE!  6, 140) ITT!  I  )  ,H( f ) , 1  =  1 , NH ) 

140  FORM  " 5( E  10.1, 2X.F10.2)  ) 

00  CONTINUE 
RETURN 
END 


(  }  ) 


onn  non  non  nnnnnn  ,-inn  oonn 


P  EL  FAS  F  2.3 


PR  EC  IP 


DA  IF 


P  31  32 


23/13/29 


SUB ROUT  I  ME  PP  EE  IP  (  NNQ,  DT  ,  A  A  ,  OR  Q,  Q3 CT  ,  A  ) 

******  ************  **  ******  **  ******** **************************  ******* 

***  THIS  SUflROUT IN  E  COMPUTES  P AI NF ALL-EXCF SS  AND  ARRANGES  THE  *** 
***  I  f !PU T  IN  A  PPDPEP  FASHION  *** 

**************************** ****************************  ************* 

D I M El  S ION  IT l(  IDO) ,  I T2(  1  30) ,01  (1  00) , TI ( 103) , K (2000) 

1 , INFI20) .  JTQ1  ( 203) ,  JTQ2I?  001  *TQ<  ?0<3)  ,008(203)  ,  INF0<  20  ) 

2.PHK103) 

COMHON/VS/TI-O) 

COMMON /VSX/X 
COMMON/VPS/PHI 

********************************  ************************************* 
***  FAF=2.54,  CONVERSION  FACTOR  FROM  INCH  to  CENTI«FTFPr,  *** 

****** ********************** ************ ********************** ******* 

FAC  =  2.  54 

****************************************************************** 

***  READ  RAINFALL  INPUT.  TIME  IS  IN  HOUP-MINUTF  SYSTEM  ,  AND  *** 

**«  INTENSITY  IN  INCHES  PEP  HOUP  *** 

***  N  NQ=  NUMBER  OF  RAINFALL  READINGS  *** 

***  S=  SOPPTTVITY,  PARAMETERS  IN  PHILIP  INFILTRATION  MODEL  *** 

****************************************************************** 

READ! 5,1)  NNQ 
I  FORMAT (15,  13X.F10.4) 

WPITFI6.il)  NNQ 

11  FOR  MAT (  5X,  '  MO  OF  RAINFALL  PE  AO INGS * , 1 5 / ) 

************************************  **************************** 


***  I NF= 
********** 

P  F  A  D  (  5  ,  3) 

3  FORMAT (20A 
WP I T  E( 6,3) 
********** 
***  PvnL= 
********** 
PFAD( 5,518 
518  FOR  M AT ( 2F  1 
P.VnL  =R  VOL* 
QVQL=0VrH* 
W  D  I  T  F  (  6,61 
618  r  no  MAT ( 5  y, 
1  ) 

********** 
***  IT1= 

******  * w  ** 

READ! 5 , 2 ) ( 
2  FORMAT ((4( 
DO  703  1=1 
700  QI ( 11=01(1 
WPITFI 6,7) 
7  FOp  MAT ( 5X , 
1 , 3X  ,  M  NT  FN 
WP I T  c  (6,8 
8  FOR  MAT ( 2X, 
1 ,3X  ,  *HnUP' 
WR  I  T  F  t  6,2) 
******  **** 
***  o  rpijc 

***  REGI 


DATE  OF  THE  EVENT  AND  THE  KATFP SHE D  IT  OCCURED  O 

************************************************ 

(  INrm  ,  1  =  1,20) 

4  ) 

( INF ( 11,1=1,20) 

****************************** ********** ******** 

VOLUME  OF  FAJNFALL  AMD  QVCL=  VOLUME  OF  SUPFACF 

********************** ************************** 

)  RVOL.OVOL 
0.4) 

FAC 

FAT 

8)  FVOL.QVOL 

RAINFALL  VOLUMF  IS  • , F 1 0. 4, 5X , • »UNCFF  VOLUME  IS 

************************************************ 

CLOrK-HOUP  ,  TT2=  MINUTF  AMD  01=  RAINFALL  INTENS 

*t** ************************************** ****** 

IT  1 (  1 ) , 1 T2 ( I ) , QI ( I ) , 1=1, NNQ) 

215, FI 0.4)  )) 

,  NNQ 
)  *F AC 

*  T  IME  •  »  3  Y»  • INTENS I TY 4X, • TIMF* , 3 X , • I  NT E NS  I T Y • ,4 
SITY* ,4X, • TIME* , 3X,* INTENSITY* /) 

•HOUR'  ,  1  X,  *MIN*  ,4*  ,  'CM /HP*  ,3X,  *HnilP*  ,1  X,  *MTN*  ,4X 
, 1 X,  'MIN' ,4X,  • CM /HR*  ,3  X, *HOUR* ,1X, *  MIN* ,4X , *CM/H 
( IT1( I) , I T2( I),QI (I ), 1=1 , NNQ) 

************** ****************** **************** 

E  CLOCK-HOUR  TIMES  TO  ABSOLUTE  TIMES  IN  SECONDS 
NNING  WITH  Z FRO 


N  *** 
****** 


*********** 

RUNOFF  *** 

********  **  * 


•  ,  F  13.4 

********* 
ITY  *** 
********* 


X, ’TIME' 


,*  CM/HR* 
R  •  /) 

***** 

*** 

**• 


C16 


(inn  non  lo.nn  non  on 


RELEASF  7.0 
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date 


P  31  32 


23/13/29 


***  CHANGF  tjmf  to  SECONDS  *** 

************************************************ ******** ******* 

Til  11  =  0.0 
DO  12  1=2, NNQ 

IFIlTim  -  IT1U-11I  14,13,14 

13  TIC  I  )=( IT21 I  1-IT2I  I— 1 1  I* 60.0 04 T1  (t -1  1 
GO  TO  15 

14  IF(ITKI)  -  IT1II-11I  IT, 17,16 

16  T  I  (  I  1=7  II  1-1  )  +  o0.0  0*  I  T2C  I  1  *160-1  T2  u-i  )  1*60.  00*360  0.30-,  IT  1C  II- 
1IT1I  I-  11-1  ) 

GO  TO  16 

17  TIC  II* TI I  1-11*60.00*  IT2I  11*3600.  00*1  T1  U  1  *60.00*160-  1T2(  1-1  i  1* 
l ( 24- IT 1( I- 11-1 1*3600.00 

15  CONTINUE 
12  CCNTINUt 

00  21  1=1, NNQ 
21  TIC  1  )=TI|1  1/60.0 
TTF IN=  T I <  NNQ 1 
HP  I  TEC  6,191 

19  FORMAT  t  7  X ,  'TINE* ,3 X,  ' I NTENSI Ty',4X,*T7ME*,3X,* I  NTENS ITY' ,4  X , *  T IMF* 
1,3X, ' I NTEN SIT Y*,4X, 'TINE ',3X,' INTENSITY* / I 

WPITEC  6,20) 

20  FORMAT  (OX,  ’MN*  ,6X,  •  CM/HR  •  ,6X,  •  MN  • ,  6X,,CM/HP.»  ,6  X  ,  •  MN‘  ,6  X  ,  •  C  M/HR  • 

1  ,6X,  'MN* ,6X,  • CM/HP '/ 1 

HP  IT  EC  6, 18 1  (TIC  II  ,QI C I) ,1=1 , NNQ ) 

19  F0PMATC4CF10.1 .F10.41I 

******  **************************  ********** 

***  COMPOTE  VOLUME  OF  PAINFALL  CCM1  *** 
********************************  ********** 

MNQ=NMQ- 1 
VP  A  I  N=  0 .0 
00  135  I =1 , MNQ 

VP  A  Jfl=  VRAIN-M  QI C  l  I*C1  1C  I  +11-TI  il  )  1  /6  0.  01 
135  CONTINUE 

HR  I T  E I  6, 140) VP  AIN 

140  FORMAT  (  5X, 'COMPUTED  VOLUME  OF  RAINFALL  C CM  1= ■ , F 1 0. 3/ / 1 

********************************************************************** 
***  P E AO  0PSEPVF0  FUNOFF  INPUT  .  time  IS  IN  HOUR-'MNUTF  SYSTE''  *** 

***  ,AN0  PiJNOFF  IN  INCHFS  PER  HOUR  *** 

***  N  QQ=  NUMREP  OF  RUNOFF  READINGS  *** 

******  ******************** ******  ************************** ********  »*** 
P  EADI 5 , 4 )  NQQ 
4  FORMAT (I  5) 

HP ITEC  6,41  )  NQQ 

41  FORMAT  (5X,  •  NO  OF  PUNI'FF  R  EADI  NG  S  *,I5/) 

********************************  ****** ********************  ********** 


***  I Nf 0=  OATF  OF  RUNOrF  EVENT  AND  THE  WATERSHED  IT  OCCUPFD  ON  »** 

******************************************************************** 

PEA0I5.5)  C  INTQC I  1,1=1,20) 

WPITEC  6,5) C  TNFQC 11,1=1,20) 

5  FORMAT  C20A4) 

*•***********•**••**••****•••******•••****•*•*•**•*••*•••* 


**♦  J TO  1=  CLOCK-HOUR,  JTQ2=M! NUTE  AND  008=  DISCHARGE  ♦** 

******  ******  ************ ************ ***•*•*•*•*•*• ******** 


READ (5, 6)  ( JTQ1C I ) ,JTQ7<  I  I .QOBCI  1, 1  =  1 ,N00) 
6  FORMAT i C4C  2I5,F10.4) 1) 

DO  705  1=1, NQQ 
705  QOBC I )=QPBC I 1*FAC 


Cl  7 


orir>  nnr\  nom  non 
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VP ITF< 6, O) 

9  FOP  M AT (5X,'T|ME',6X,  'RUNOFF* ,4 X,  'TIME*  ,  6  X  ,  'RUNOFF* ,«», *T IMF* ,6X, 
I'PUNlirF'  ,4X, 'TIME*  , 6 X,  'RUNOFF  •  /) 

WP I TE( 6, P ) 

WRITE! 6,6)  ( JTQ1 (  I  )  , JTQ2 I  I  ),QOBI  I)  ,I=1,NQQ) 

******************************** *» 

***  CHANGE  TIME  TO  SECONDS  * ** 

**********  **********************  «* 

TQI l )=0.00 
DP  55  1=2, NOO 

IF(JTC!(I)  -  JTQKI-1))  52,51,52 

51  TCI  I  )=  !  JTC2II  )-JT02< 1-1) J*60+TQI  I- 1 ) 

GO  TO  53 

52  IFIJTCim  -  JT01(I-m  56.56,54 

54  T0(  I  )=  TQ 1 1  - 1  I+60PJTQ2I  I  )  +  (  60-J  TC2I  I-l )  )  *60  +  3600*  (  J  TO  1 !  II-JTQKI- 
ll)-l) 

GO  TP  53 

56  TO!  I )  =  TQ(  I  -  1 !+60*JTQ2(  I ) +3600* JT Ql II  1  +  60*160— JT 02 1  I-l) 1  +  I24-JT01 
1(1-  11-11*3  600 

53  CONT INUF 

55  CONT INUF 

***************************************************** 

***  SPPT  OUT  THE  PEAK  RATE  OF  RUNOFF  OBSERVED  *** 

***  P  SO=P  E  AK  RUNOFF  *** 

********************************  ********************* 


JP=1 

OBO=OOB( 1) 
on  25  1=2, NQO 
TO!  I )=  TO!  I  ) / 6 3 . 0 
IFIPPQ.GT.OOBI  in  GO  TO  25 
ORQ  =  QOB(  I) 

Jp=! 


25  CONTTNUC 
OBQT  =  T  0 ( JP  ) 

WRITE! 6, 27) 

2  7  FPr  MAT ( 7X,  •  T IMF  •  ,  5  X,  »P  UNOFF  •  , 3 X ,  'TIME* ,5X,' RUNOFF*  ,6X 
ISUMPFF ' »  6X  « ' T I MF  ' « 3X,'RU"PFF ' / ) 

WR I T FI  6,26) 

26  EPPMAT ( 8X,  ' *  , 6X , 'Cm/HR • , 7X , • MN • , 6X , • CM/HP •  , 7 X , • MN* , 
IX, '"N* ,6X, 'CM/HR'/) 

WR ITT! 6, 18  1 1  TO ( I), 00 R( I ) ,1 =1 , NQQ ) 


, »T !«E'  , 6  X 


6  X  ,  'FI*/  HR  • 


•**  CnfPUTF  VOLUME  OF  SURFACE  RUNOFF (CM)  *** 

****** ******** ************************ ********* 


VRUNOF  =0.0 
on  145  1=2, NOQ 

VPUMOF  =  VFUNOF  +  (QaeiI l+QOPI I —  1 )  )*0.  5* (TQI IJ-TQI I-l) 1/60.0 
145  CONTINUE 

QVOL=V  RUNO  F 
WRITE!  6,  150)  VPUMOF 

150  rnpMAT ( 5X, 'COMPUTED  VOLUME  OF  SUPFArE  RUNOFF (C M )=• ,F 1 0 . 4/ ) 

************************************************** 


***  change  punofr  to  cubic  metfrs  pfp  second  *** 

********************************  ********  ********** 


DO  715  1=1, NOO 

715  QnR( I )=QPH( I )*A*100. 0/36.0 
WRITE! 6, 27) 

VP ITEI 6, 720) 


7 
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720 


115 


125 


69 


130 

120 

627 

no 


91 

93 


92 


1  ™P'.'1£I,<.824  '  C  UM  /  S  •  ,7X.*MN« ,6X,*CUM/S*  ,7X,  *MN*  ,6X.*CUM/S', 

i  *  *  t  *  i  oX  •  *  C IJf*  /  S  *  /  l 

WRITE!  ft,  lfl  )  f  T  Of  IJ  ,QOB! I  1  ,  1  =  1,  NO  01 
**************************************************************** 
C0MPUT E  INFILTRATION  LOSS  .AND  SUBTRACT  IT  FPOM  RAINFALL 
************* * ************ ****** **  ** ******  ********************** 
00  115  I=1,NN0 
Tim  =  TI(li/60.0 

3Ne?S?ES!??8 

CALL  INFIL (NMO.S.AAl 

******  ************  ********************  ************* 

***  DETERMINE  THE  PAINFALL-EXCE  SS  VOLUMEICM)  *** 
******************************** ** ***************** 

VFXCES=C 
00  69  1=1, MNQ 

P.LOSS=  (PHI  (  I  1*RHI(  I*  111*0.  5 
Q I  (  I )  =  QI ( I J-PLOSS 

ifiokii.lt. o.oo)  oim=o.oo 
temp=t  it  i*n-Tn  i) 

VEXCFS=VEXCES+QI ( I )*TEMP 

CONT INUF 

Qll  NNQ 1=0.00 

NNQ=NNQ-H 

T I  (  NNO 1  =  T I  I NNQ-1 )  ♦  0.  5 
01 1 NNQ) =0.0 
WRITFI6, 130)VEXCES 

’COMPUTED  VOLUME  OF  RAINFALL  F  XC  ESS  (  CM)  =  •  ,  FI  0 .4/  ) 

OP  1 20  ;=lfPNO 
Til  I )  =  T 1 1  1  1  *60.0 
WRITE!  6,627) 

Ff'PMAT(5X,  'RAIN  INFILTRATION  IS  AS  FOLLOWS  *//) 

WRITE!  6,110) 

iMNFtLT!is;;i!^EC5x:J!Nn:!/),,TiME,’4x,,,NPiL,,6)<’,TiM-,’5x’ 

WRITE!  6,26) 

WP  ITC!  6,  18)  I  TIC  11  ,  PHI  U)  ,  1  =  1,  NNQ) 
************=***•***************  * 

***  ADJUST  THF  TIME  SCALP  *** 
********************************* 

IF( 0  I !  1 )  -  0.00)  92,92,91 

IF!  TM  NNQ )  -  TTFIN)  93,93,90 

NNQ=NNQ»1 

T I « NNQ ) =TTFIN* 10. 00 

QI  !  NNO ) =0. 00 

GO  Tri  90 

CONT  INUC 

IK=0 


I FR  P  =0 

DO  76  1=1, NNQ 
IN*  I 

PTMP=QI<  II 

IF  I PTM  P-0. 00 )  76,76,  77 
77  I K= I K*  1 
TMP  =  T I  (  I  ) 

GO  TO  78 


♦  *** 
*** 
**  ** 
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76  CONTINUE 

78  CONTIN'JF 

I F  (  IK.CC.O  )  I  FP.P  =-  1 
IF  I  IEPP.ro. -1  )  GO  TC  23 
DC  79  I = I N ,  NNO 
PTMP=Q I ( I ) 

Oil  IK) =  PTMP 
T  1  *  IK»=TI( I )-TMP 
IK  = IK*  1 

79  CONTINUE 

OU  I K  )  =0.0  0 

Til  1 K )  =  T I (  !K-1 H 10-37 
IK* IK+  l 

on  ik)  =0.00 

Til  IK) =T I (  I K— l)vl0.00 
NNQ=  IK 
90  continuf 

******  *************************************** 

***  SHUTOUT  THE  PEAK  PATE  OF  PAINFALL  *•* 

***  Q IMAX  =  PEAK  PAINFALL  INTENSITY  *** 

********************************************* 

QlMAX=OH  1  ) 

no  22  1=2. NNO 

IF! HIMAX.GT.Oiri >!  GO  TO  22 
QIMAX=QI < I ) 

22  CONTINUE 

HP  TTEC 6. 26 )  OIMAX.OBQ 

24  FORMAT ( 5X,*PFAK  RAINFALL  INTENSITY  I  CM/HP) *, F10.4, 2X ,*PEAK  RUNOFF 
1 ( CM/ HP )*,F 10.4 ) 

23  CONTINUE 

HP  I T  E(  6,  62  8  ) 

628  FORMAT < 5X, 'PAINT AL L-EXCESS  IS  AS  FOLLOWS  ' // ) 

WRIT!  6,  19) 

WPITFI 6,20) 

WPITFI  6,1" )(TI (I QI  Cl ), 1=1, NNO) 

WPI TT( ft, 160)  TUP 

160  FOP  MAT ( 5X ,  'LAG  T IMF  BETWEEN  THE  START  OF  PAINFALL  AND  T  HF 
1FFFFCTIVF  PAINFALL  < M I N ) = • ,F l 0 .2 / 1 
RTTURN 
FND 
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SUBROUTINE  XDATA<NNQ,NXM,DT) 

***************************************************************** 

***  THIS  SUBROUTINE  COMPUTES  THE  RAINFALL  EXCESS  AT  EQUAL  *** 
***  T I  ME  INTERVAL  *** 

****************************** ** ********************************* 
DIMFNS ION  Cl  I  100  I, Til  1001 , XI  20  00  ) 

COMMON / VS/ T I, 01 
COMMON /VSX/X 

************************  ********  ******* 

***  DEFINE  LIMIT  OF  COMPUTATION  *** 
********************************  ******* 

DO  40  1=1, UNO 
TH  I>=  Tit  IJ/60.0 
40  CONTINUE 
NMAX=NXM 
QIINNQ *11=0.0 
TI!NN0*1)=TI(NNQH500.  0*DT 
NXM=T I (NNQ l/DT+l.O 
IF( NXM.LT .NMAX 1  NXM=NMAX*5 
TC=0.0 
KK=  1 

00  60  1= 1»  NXM 
IFITC.GT.TKKKIl  GO  TO  65 
GO  TO  TO 

65  xm=C  (TI(KK  l-(TC-OTn*CI(KK)*<TC-Tr  (KK»I*QI  (KK+IU/OT 
KK=KK* 1 
GO  TO  75 
70  Xt I  1=0  riMK  1 
75  CONTTNUC 
TC=TC*OT 
60  CONTINUE 
RETURN 
END 


C21 
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SURPOUT IMP  NEWTONINNQ, YEXC ES,X3,R) 

*********************************************************** 

***  THIS  SUBROUTINE  CnMDUTES  THE  SOPPTIVITY  RAPAMFTEPS  IN  THE  *** 
***  PHILIP  TWO-TERM  I  NF  1 1  T°  A  TI ON  MODFL.  *** 

***£*********#**********«***«***  ***«**•**»•****«*• ********** ********** 

DIMENSION  TKIDOI.QI  (lOOl.PHK  1001 
COMMON /VPS /PHI 
COMMCN/VS/TI,QI 
X 1=  0 . 1 3 

X  2=  2 . 5 
K0Uf!T=3 

10  CONTINUE 

CALL  I NFIL (NNQ,X1 , B) 

MNQ=NN 0-1 
VEXCES  =0.00 
00  60  1=1, MNQ 
AVG= ( P  H I ( l  I*PH1|1  +  1)1*0.5 
OIFO=0  im-AVG 
IFIOIFQ.lt. 0.001  DTFQ=0.00 
VEXCES  =  VFXCES«-D1FQ*< TI( I ♦!  l-TIIl  11 
60  CONTINUF 

FX  1  =  VE  XC. ES-YEXCE  S 

CALL  INFTL INNQ,X2,B) 

mnq=nno-  1 

VFXCES=0 .00 

DO  55  1=1, MNQ 

AVG  =  ( P HI  (  I  ItPHK  1+  11  1*0.5 

OIFO=Qtm-AVG 

IF! DIFQ.LT .0.03)  DIFQ=0.00 

vExcEs=vrxcES+oiro*(  Tin  4-1  1-ti  nil 
5*  CONTINUE 

FX2=VEXCFS-YEXCES 
OF X  =  FX  2-rX  1 

IFIDFX  .NT. 0.031  GO  TO  6 
X3=  2 .0  *X  2 
GO  TO  7 

6  CONTTSMjr 

X3=< Xl*f X2-X2*FX1 1/OFX 

7  CONTINUE 

ppn  =  AB  S 1 IYFXF'  S-VFXCrSl/YEXCES) 

KOUNT  =  KOUNT*  t 

1F{ PRP  .LT. 0.01 1  GO  TO  25 

IF  (FOUNT. GT. 201  GO  TO  30 

XI  =  X? 

X2=X3 

GO  TO  10 
30  CONTINUE 
KPITFl 6,36) 

36  FORMAT (5X, 'ERROR  IN  THE  LOOP  AND  NO  CONVERGENCE' /) 

25  CONTINUF 
RETURN 
END 
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SUBROUTINE  I MF I L  C  NNQ  *  S  «  A l 

********************************************************************* 

***  THIS  SUBROUTINE  COMPUTES  INFILTRATION  USING  THE  TWO-  TERM  ♦** 
***  INFILTRATION  MODEL  *** 

**********  ********************** **  ************ ************ ******  ***** 

01  MENS  I  ON  THIOOItPHIC  1001 
COMMON/VS/T! 

COMMON/VPS/PHI 

******  ******************************** 

***  TWO  PAR AMF TF PS  ARE  S  AND  A  *** 

******  ************************  ******** 

00  5  I  =  2» N NO 
TEMP  =  T  1(1) 

PHI  (1)  =0.5*S*TEMP**( -0.5)  * A 
5  CONTINUE 

PHI  (  1)  =PHI  ( 2) 

RETURN 

END 


C23 
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suppojt  ir:r  pm  sen  (  n.mst.ipt.fps.nobs) 

DIMENS  I  Of'  E(  10)  ,B(  10)  .C(  1  o),D(  10  )  ,  IS  (1  0)  ,X  (10)  ,YI  10)  ,Z  (  10)  , 

1XM(  10)  ,XL(  10), S( 10,101 ,HP I  10) 

CnMMON /RA°A/X,XL,XM 
COMMON /  Pr>SE/Y 

•  ******:************♦******#********** 


THIS  PROGPAM  SOLVES  CONSTRAINED  MINIMIZATION  PRORLEM,  WHICH  *** 

REQUIRES  THAT  THE  VECTOP  ALWAYS  RE  AN  INTERIOR  POINT  OF  THE  *** 

EFASIPLE  SET  *%* 

THE  SOLUTION  TECHNIQUE  IS  A  MIXED  APPLICATION  OF  ORIGINAL  *** 

ROSENBROCK  METHOD,  PALMER  VERSION  AND  PENALTY  FUNCTION  METHOD  *** 

THE  USER  MUST  SUPPLY  A  SUBROUTINE  ORJFCT (VALUF • NOBS, N)  FOR  *** 

EVALUATION  OE  THE  OBJECTIVE  FUNCTION  *»* 

THE  USER  MUST  ALSO  PROVIDE  THE  FOLLOWING  INFORMATION  *** 

N=  NU*REP  OF  VARIABLES  *** 

X(I)=  INITIAL  GUESS  OF  THE  VECTOR  *** 

XM(I)=  UPPER  LIMIT  OF  THE  PARAMETER  XIII  *** 

XLII)=  LOWER  LIMIT  Of  THE  PARAMETER  XIII  *** 

EPS  =  CONVERGENCE  TOLERANCE  BASED  ON  THE  CHANGE  OF  OBJECTIVE  *** 

FUNCT  ION  *** 

MST  =  MAXIMUM  LIMIT  OF  NUMBER  OE  STAGE  SEARCH  *** 

I PT  =  3 - ONLY  the  FINAL  ANSWEP  IS  DESIPEC  TO  BE  PRINTED  *** 

TPT=  1  -  INTERMEDIATE  VALUES  OF  EACH  STAGE  SEARCH  IS  DESIRED  * ** 

TO  BE  PP I NT  ED  *** 

I PT=  2 - INTERMEDIATE  VALUrS  OF  EACH  CYCLE  SEAPCH  IS  DFS I  RED  *** 

TO  BE  PPTNTFP  *** 

SET  INITIAL  STEP  LENGTH  ANO  DIRECTION  *** 


ft******************************************************************* 


DO  1  I  =  1 ,  N 
0(1  )=0.)*X( I ) 

Y(  I  )  =  X  (  I  ) 

00  2  J  =  1 ,  N 
5(1  ,  J)=0. 

IF(J.EQ.I)  S(I,J»=l. 

CONTINUE 

CONT INHE 


STAFTING  Of  STAGE  SEARCH  *** 
*************-*  ******  **«******« 


NEE  =  1 

CALL  OBJFCT(VALUE,NOPS,N! 

P0= V  AL  UE 

3  NS  =  N  S*  1 
OBJ=RO 

DO  4  1 =1,N 
NR(  I  1  =  0 
I S (  I  1=0 
Z( I  1=0. 

4  CONTINUE 

IF< I PT .EO. 0)  GO  TO  5 
WP I T  E( 6,650) 

650  FOP  MAT ( /40X, 40H***** ******************************** »*» ) 

WPITF(  6,604)  NS 

604  FORMAT ( //4BX, 1 8HSTAGE  SEAPCH  - ,15) 

*******»«*******•****•*••»*«*•««  ********  ****************  ******** 

STARTING  OF  CYCLE  SEARCH  (WITHOUT  CHANGING  SEARCH  DI R  FCT ION)  *** 

****************************************************************** 
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5  NC= - 1 

6  NC=NC*l 

DO  7  1=1,  “1 
I C  0  fJ  =0 
ox*  n  1 1 > 

32  on  8  j  =  i.n 

yi  j)=x  i  j)«-nx*s(  i  ,j> 

********************  ******** 

CHECK  IF  IT  IS  AN  INTFRIOP  POINT  "F  THE  FEASIRLE  SET  *** 

********************************************************** 
IF(  YIJ).LT.XL(J)  .OP,.  YIJI.GT.XHIJII  IC0N=1 
8  CONTINUE 

IF!  ICON. HE .  1  )  GO  TO  25 
0MIN=1 .OHIO 
DO  26  J=1,N 

IFISII. JI.FO.D.)  GO  TO  ?6 
IFIYIJl.GT . XL ( J I )  CO  TO  30 
DD= ABS I  I  XI Jl-XLl J11/SII,  J) I 
T  F  I  DD.LT.DMjri)  DMIN'DD 
30  IFIYIJI.IT.XMIJH  GO  TO  26 
DD= AOS ( ( XI JI-XMI J))/S( I , J»J 
IFIDD.LT .DMINI  DM! N=DD 
26  CONTINUF 

IFIOX.LT.O.)  DM  I N=-0MI N 
DX=  DMI N 
ICON=2 
GO  TO  32 
25  NEF=NFF*1 

CALL  OBJECT! VALUE, NOBS, Ml 

********************************************************** 

CHECK  IF  THE  VALUE  OF  OBJECTIVE  FUNCTION  IS  INOPOVED  *** 

********************************  ************************** 

IF(  VALUE. GT.t>UI  GO  TO  9 

************ ******••*•«* **************************** •***•*• 
THE  VALUE  OF  OBJECTIVE  FUNCTION  IS  IMPPOVED  (SUCCESS)  *** 

****************************** **•***•**•****•••* *********** 


pn=  VALUE 

I F (  ICCN.ro. 2)  GO  TO  37 

IF  (  IS<  I  1.NE.2)  IS!  I  )  =1 
Z(  I  1=2  I  I  )*D(  I  1 
0< I  1  =  3  .*01  I  ) 

C,r  TP  31 

33  Zt I  )=Z II) »0V!N 

0(  I  )=-DMIM-0.5*Dt I  ) 

IS(  I  1=2 

31  DO  10  J  =  1 , N 
XI  J  )=v  (  J  1 
10  CONTINUE 
GC  TO  7 

******  ************ ************  ******** ******•*•*•• **•*•**••*** 


THC  VALUE  OF  OBJECTIVE  FUNCTION  IS  NOT  I MPRDVFD  (FAILURE)  *** 

*•*****••**•***«**••*****••••*******•***••******•*•**••*••**•* 


9 

7 


IFIlSin.EO.n  I  SI  1 1  =2 
DI I  )=-0. 5*01  T  I 
CONT INUF 

************** ••*•**•**** 

END  or  CYCLE  SEARCH  *** 

************** *********** 


C25 


onno  ono^  om  ^ono 


PELEASE  2.0  8P0SEN  DATE  *  <»31?2  23/13/2° 

IF(  IPT  .ME. 2)  GO  TO  11 
WPITEI  6,  605)  NC 

635  FORMAT  I //48X, 1RHCYCLE  SEA°CH  - ,15) 

WP ) TF(  6, 60  6  )  P0 

606  FORMAT ( //34X.32HTHE  CURPFNT  OBJECTIVE  FUNCTION  =  ,  E  20 . 8  /  /  60  X ,  2 1  HT  HF 
1  CUPPEMT  VECTOR  IS) 

WP I T  Et  6,60?)  (XII), 1=1, N) 

******  ******  ********* 

CHFCK  IF  ALL  the  SE4PCH  OTPECTION  HAVE  FAILFD  AT  LEAST  CNCF  *** 
FOLLOWED  BY  A  SUCCESS  *** 

******  ************** ****** ******  ****************** ******** ******* 

11  I  FT  =  0 
on  12  T  =  I ,  N 
IFl  I  S(  I  )  .ME. 2)  IPT=1 
IF( ISC  It. EC. 01  NR  C I ) =NR ( 1 1 ♦ 1 

12  CONTINUE 
DO  34  1=1, N 

IFCNPCD.LT.A)  GO  TO  35 
24  CONTINUE 
Gn  TP  14 

35  IF (  IPT .FO.  1)  GO  TO  6 

************************* 

END  OF  STAGE  SEARCH  *** 

************************* 

IF  I  IPT .FO. 0)  GO  TO  13 
WP  I  T  El  6, 65  3 ) 

NC=NC*  1 

WP  I  T FI 6,607)  NC 

637  FOPMAf ( //41X.30HTOTAL  NUMBER  OF  CYCLE  SF  ARCH  =,!5) 

WP.  I  T  EC  6 , 60  6 )  PH 
WRITF(6,603)  (X(I),I=1,M) 

*********************************************************************** 
CHECK  IF  THE  RESULT  IS  SATISFIED  WITH  THE  PPEASSIGf,cD  CONVERGENCE  *** 
TOLFRANCF  *** 

********* r************************ **************** **********  *********** 

13  ircconj-pm.Lr.EPS)  go  to  14 

********************************************************************* 
Check  it  THE  NU'-’BCF  OF  STAGE  SFAPCH  GofiTpR  than  ASSIGNED  LJMIT  *** 
****************************************************************** *** 
1FCNS.LT.MST)  GO  Tf)  15 
WPITEI 6,650) 

WPITEI 6,608)  MST 

6  OR  F0PMAT(//40X, 18HDO  NOT  CONVERGE  IN , I  5, 5X , 1 4H ST  AGE  SEARCHES) 

WPITEI 6,606)  PO 
WRITFI 6,6^3)  IX< I),I=l,N) 

603  FORMAT l2ri0.4,Fl0. I) 

GO  TP  700 

********************************************************** 

CALCULATE  NEW  SEARCH  DIRECTION  FOR  NEXT  STAGE  SEARCH  »** 

PAL  mtp  S  VERSION  IS  USEP  TO  COMPUTE  THE  NEW  DIPECTION  *** 
********************************************************** 

15  DO  17  1  =  1,  N 
SUM  A=0  . 

DO  18  J=  1 , N 
El  J ) =0  . 

DO  IQ  K=  I  ,  N 
F(J )=F (J )+Z(KI*S(K,J I 
IP  CENT  INUF 


C26 


ono 


P  EL  FAS  F  2.0 


RPOSEM 


04  TC 


8313? 


?  3/  1  3/29 


SUM  A  =  S  UM  A*  E  (  J  )  2 

18  CONTINUE 

AA=S0PT|  SUMA) 

IF( AA. EO.O  .  )  GO  TO  3 
l F (  T  .EQ.  1)  C.D  TO  20 
IFlAOSm  I-l)  I.LE.FPS)  GO  TO  27 
0A=  1  ./ SORT l AR**2-AA**2) 

RA= AP/ AA 
CA=0A*R  A 
CB=D4/RA 
DO  21  J=1,N 
C(  J  )=S  (  I  ,  J  ) 

S(I.J1*F(J  )*CA-B( J )*CP 

rijUfui 

21  CONTINUE 
GO  TO  23 

20  no  22  4=1,  N 
C(J)=S<!,J) 

S(ItJ)=F(J)/4A 
BI  J  )=F  (  J  ) 

22  CONTINUE 
GO  TO  23 

27  DO  23  J*1,N 
CTC  M=S t  I  ,  J  ) 

S(  1  •  J  I  =C  ( J  ) 

C< J  >=CTFM 
BI  J  )=E  l  J  ) 

23  CONTINUE 

23  AB=  A  A 

17  CONTINUE 
GO  TO  3 


A  MINT  MUM  HAS  REEN  FOUNO  *** 


14  NS=NS*l 

won  t.;  8,630) 

H°  IT  M  6,  60  9)  NS,  NCF 

609  FrO'*ATI//4PX,24HA  MINIMUM  HAS  BEEN  c OUND  / /4 1  X ,  3 9HT  OT  AL  NUMRFP  OF  S 
1TAGF  S  CATCH  =, 13//  39X.37HTOTAL  NUMBER  OF  FUNCTION  EVALUATION  =,15) 

WPITr( 6,610)  PO 

610  FOP.MAT  I //3PX,  23HOPTI  MI  2AT10N  FUNCTION  =,E2  0.  8//50X«15HFI  NAL  VECTOR 

WPITTi!  6,6031  (XII), 1  =  1, N) 

700  C0NTiNUE 
RETURN 
END 
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noon  noono 


RELEASE  2.0 


c ON VOL 


DATF 


R  3  1  3  2 


23/ 13/23 


SUPRruT  INF  CnNVOLI  PT ,NX«, 0P,QPT,  KOPT ,A) 

•**  THIS  SUBROUTINE  C.DNVOLUTFS  T HF  PAINFAIL  FXFFSS  WITH  *** 

***  THF  tnst antaneuus  unit  hydpogpaph  TO  DETERMINE  the  *** 

***  DIRECT  RUNOFF  HYDP  OGRAPH  *** 

•  ********************************* 

01  MENS  IC.N  HI  2000)  ,  X<  2000)  ,  Y I  2000  )  ,  TY  I ?  00 0 ) 

COMMON /T  HP FF/H.NH 
COMMPN/VSX/X 
KNU  MP=  2000 
NHA=NH+ 1 
NMP  =  NH +f'y  M 

IFINHB.GT.KNUMB)  nhb=knumr 

NYMsMHB-X 

00  30  l=NH A , NH  R 

hi  n=o.p 
30  continue 

MXM=NX  M+ 1 
on  10  I=MX  M*  NHB 
XI  I  )=0 .0 
10  CONTINUE 

tc=o.o 

VMX=0. 0 

Y  B  =  0 . 0 

on  3  5  N= 1 ,  NYM 
tc=ti>ot 

YIN)=0.0 
TYIN)=TC 
J  =  N 

IFIN.GT.NXMl  J=NXM 
00  40  1  =  1,  J 
I  A=N  +  1-T 

Y ( N  )  =Y I N  )  ♦X  I  I  )*Hl I  A) 

40  CONTINUE 

Y(M  )=Y IN 1*0T 

IFIYNX.lt.  YIN'I  |  YMX= Y  IN) 

YMN=0.  0 1  * Y  MX 
NOY=N 

IFIN.GF. KNU MR )  GO  TO  45 
IFIYIN ).LT .Y  MU . A  NO .N  «GT. 20 )  GO  TO  45 
35  CONTINUE 
45  CONTINUE 

00  5  1=1, NOY 
TFMP=TYI NO Y* 1- I ) 

Tv ( NOY-l+2)=TEMP 
5  CONTINUE 
TYl l)=0.0 
NOY  =  NO  Y  +  1 
DO  15  1=1, NOY 
15  TYI  I  1=  TY I  I  ) *60. 0 

******************************************************************* 

***  SOFT  OUT  THE  PEAK  RATE  OF  SURFACE  PUNOFF  AND  ITS  TIME  OF  *** 
***  OCCURRENCE  *** 

******************** ************  ****************************** ** *** 

QP= YU) 

On  53  1=2, NOY 
iriOP.GT.YI I))  GO  TO  50 
0P  =  v ( I  ) 
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ouu 


RELEASE  2.0 


CONVOL 


niTr 


B312Z 


23/ 13/?Q 


!R=  1 

■SO  CONTINUE 
QPT  =  TY ( ID) 

IF  (  K  OP  T  .  NF  .  0)  GO  Tfl  60 

WRITE! 6, 145) 

145  FOP R AT ( 53,  'THE  SURFACE  RUNOFF  HY OR OOP A PH  IS  AS  FCLIOWS'//) 

WP  ITF! 6,  150) 

150  FORMAT ( 7X,  • TJ YE* ,6 X, 'RUNOFF*  ,5X,  'TIME* ,5 X, 'RUNOFF'  ,7X,'T1ME',5X,' 
1  P.UNHFF  • , 5X ,  'TIME ',  5X,  •  RUNOFF  •  /  1 
WRITE! 6, 155) 

155  FOPMAT l 8X, 'MN» , 7X, 'CM/HP • ,7X,'MM ', 8X,'CY/HP' ,7X, 'YN* ,RX, 'CY/HR' , 7 
IX,' YN'  , 8  X ,  'CY/HR'/) 

WR  I  TE(  6,  16  0  )  (TY(J)fY(J)  ,J=1,‘)0Y) 

160  FCPMAT(4( F  10. I . 2X.F10. 3)  I 


*** 


CHANGE  RUNOFF  TO  CUBIC  YF  TF  R  S  DFR  SFCOND 


FAC=A* 100. 0/36.0 
00  170  1=1, NOY 
170  Yfl  )=Y  !  1  )  *  F  A  C 
WRITE! 6, 150) 

WRITFI 6, 175) 

175  FORMAT ( RX,  'MM'  ,7X, 'CUM/S' , 7X, «  HN  « , 8X ,'CUM/S’  ,7X, 'UN'  ,PX, 'CUY/S* 
l,7X,'MN*,7X,  'CIJM/S'/) 

WRITE!  6, 160)  I TY! J ) , Y<  J)  , J  =  1 ,NOY) 

60  C0NT1NUF 
RETUPN 
ENO 
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